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1.0 INTRODUCTION AND SUMMARY 


The Orbiter Ku-band radar is presently being developed [2] 
the hope is that the present design is "molded in concrete" as much as 
possible. With that in mind, when problem areas arise in the radar 
performance, the solutions recommended in this report are those which 
involve minimum cost, minimum (essentially zero) design changes in 
hardware, and little variation in software development. 

The performance of the Ku-band radar is analyzed in detail in 
[1], and the performance is updated and summarized in sections 2.0 and 
3.0 of this report. In so doing, two different radar design philosophies 
are described, and the corresponding differences in losses are enumer- 
ated. The resulting design margins are determined for both design 
philosophies and for both the designated and nondesignated range modes 
of operation. In some cases, the design margin is about zero, and in 
other cases it is significantly less than zero (-4 to -5 dB). With 
the point of view described above, the recommended solution is to allow 
more scan time (2 minutes as opposed to 1 minute) but at the present 
scan rate. With no other changes in the present configuration, the 
radar will meet design detection specifications for all design philo- 
sophies at a range of 11.3 nautical miles. 

In the tracking mode, the Ku-band radar meets specifications 
satisfactorily at all ranges in the range and range rate tracking 
loops. This is with the assumption, however, that there is no signi- 
ficant degradation due to angle tracking errors. As is shown in 
sections 3.0 and 4.0, this is not the case. At short ranges, the angle 
tracking loop and angle rate estimator are not reliable with the present 
configuration The reason for this is primarily due to target glint 
effects which can and do become significant when the target fills more 
than 10% of the 3 dB beamwidth of the antenna pattern. In addition, 
the monopulse difference patterns are unpredictable at short ranges. 

In the next section, the recommended solution to this problem is 
described. This solution is effective in all dimensions of the point 
of view described above. 
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1.1 Recommendation for Angle and Angle Rate Tracking 

Problem at Close Ranges ~ 

A rough description of the suggested solution to the proximity 
operations problem in the angle tracking loop of the Ku-band radar is 
presented. The approach centers around slaving the Ku-band antenna to 
the TV camera during a portion of the approach phase in the rendezvous 
and docking of a target. 

The situation with the present configuration is as follows: 

There is a transition period during approach in which neither the 
Ku-band radar nor the two TV systems individually provide all the 
desired information at acceptable accuracies to carry out a rendezvous 
and docking profile. Using the two TV systems anc[ the Ku-band radar, 
however, all of the necessary position and velocity information can be 
satisfactorily obtained. 

There is presently a TV camera at each end of the Shuttle's 
stowage bay. Each camera has an 8:1 zoom lens and is mounted on a pan 
and tilt, thereby giving each camera two degrees of angular freedom. 

In addition, each pan and each tilt is encoded into a video signal 
using a Binary Gray Code of 11 bits. This corresponds to angular inter- 
vals of 0.18 degrees. The accuracy off each shaft has been quoted as 
0.3 degrees. These four encoded angular signals are readily available 
at the pans and tilts. The operations interface (OI-MDM) does not 
presently receive this angular information, and the information is not 
now being telemetered back to a ground station. However, these video 
signals with the angular information can easily be sent to the OI-MDM 
interface. This is shown as connecting lines (a) in Figure 1. Also, 
both TV cameras can presently be simultaneously monitored by the 
rendezvous operator on two adjacent CRTs. 

The Ku-band presently "talks" to an MDM, as shown in Figure 1. 

As described in sections 3.0 and 4.0, the angle and angle rate accura- 
cies of the radar become unacceptable at a larger range than do the 
accuracies of range and range rate. When the radar provides no useful 
angle and angle rate information, however, the antenna must clearly 
be kept pointed toward the target in order to provide accurate range 
and range rate information. 




Figure 1. General Block Diagram for Slaving Ku-Band Antenna to Two TV Cameras at Short Ranges 
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The suggested solution would slave the Ku-band antenna to the 
TV cameras, which are manually pointed toward the target. Both TV 
cameras would be used, and a relatively simple software algorithm 
would be developed to best use this information to direct the Ku-band 
antenna toward the target. To accomplish this, a communication link 
(cable) would have to be created so that the MDM and OI-MDM can "talk" 
to each other. From a hardware implementation and cost point of view, 
this appears to be the largest item and should be nominal compared to 
all the other alternatives that have been examined (primarily the addi- 
tion of a laser or lasers). These alternatives are described in section 
5.0 of this report. 

In summary, the suggested solution requires cable and software 
development. Neither of these, however, appears to be extensive. 

1.2 Ku-Band Radar Simulation 

A model of the Ku-band radar has been developed and is presented 
in section 4.0 of this report The model is described such that it is 
directly programmable into a digital simulation. It can operate in "real 
time" with a man in the loop and with arbitrary maneuvers in the orbital 
geometry. Some minimal computation rate will be necessary, however, in 
order to operate in real time. 

The model is given as a simulation modular package so that it 
can be placed in any of the present Orbiter simulations at NASA JSC 
with minor interfacing. 

One of the purposes of the digital simulation of the Orbiter 
Ku-band radar is to determine operating procedures when operating manu- 
ally. This is particularly pressing with the present radar schedule 
due to the Skylab trajectory £13]. 

The successful completion and operation of the proposed radar 
simulation should add significant confidence to the concept that satel- 
lites can indeed be tracked and docked and a complete rendezvous per- 
formed. Some of the example performance computations are carried out 
for the LDEF Satellite. 
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2.0 TARGET ACQUISITION - PASSIVE TARGET 

A significant part of the overall Droximity operation of the 
Ku-band radar is the capability to detect and re-acquire the target at 
all ranges less than 10 nmi. The acquisition capability of the present 
configuration of the Ku-band radar is summarized in this section. The 
development in Appendix A of [1] is used extensively. 

There are several facts which are used as assumptions and taken 
advantage of when assessing the acquisition capability of the Ku-band 
radar. These assumptions are actually assets, in particular: 

(a) The target size is known . Other properties are known, such 
as shape, composition, and possibly even the average radar cross-section 
(RCS) This provides preliminary information about expected target 
strength at the time of acquisition. 

(b) The target's angular rotation is small . This is a self- 
evident requirement if a rendezvous and eventual placement into the 
Shuttle payload bay is to be carried out. 

(c) There exists no jamming . The passive targets will be truly 
passive, and the active targets will be cooperative. No additional 
signals of an interfering or jamming nature are anticipated. 

(d) In almost all configurations, there exists no clutter . The 
only possible exception is when the target is directly below the Shuttle 
and the radar sees earth clutter. If clutter does present a problem, 
various methods have been suggested to overcome it, including the appro- 
priate rate of frequency hopping. 

(e) One possible situation is AGG saturation in the detection 
mode at close ranges. This is expected to be highly unlikely because 
detection at close ranges invariably will be a reacquisition. In that 
case, the TWT is either attenuated or removed so that the power level 
will not cause AGC saturation during at least a portion of the short 
range detection time. 

The primary emphasis in the evaluation of the Ku-band radar during 
acquisition is "how the radar performs." Comparison will be made to 
the specification values, but no heavy emphasis is placed on it. In 
those situations where the performance does not meet the specification 
requirements, reasonable alternatives are proposed which will provide 
satisfactory performance in all rendezvous profiles. 
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2.1 Target Detection Processing Losses 

The target detection processing losses as they are presently 
estimated are summarized in Table 1. The losses anticipated by both 
Hughes and Axiomatix are shown. The reasons for the agreements and 
disagreements and the resulting philosophies in radar performance are 
as follows: 

The transmit loss is a hardware loss. It is presently estimated 
at 3.7 dB by Hughes, and Axiomatix believes this is a realistic value. 

The beam shape loss is listed as an average of 2 dB by Hughes 
and 3.2 dB maximum by Axiomatix. The difference is seen by inspection 
of Figure 2. If a target remained at boresight during its entire period 
of illumination, the antenna beam shape loss would be zero. The target 
location for maximum beam shape loss is also shown in Figure 2. The 
loss for this path is the one used by Axiomatix. If the beam shape 
loss is averaged over the region indicated m Figure 2, the result is 
approximately 2 dB when the scan overlap is 15% to 30%. This average 
is reasonably insensitive to the amount of overlap when the scan overlap 
is within this region. For a detailed analytical description of beam 
shape loss, see [4]. 

The pre-sum mismatch has been evaluated as 0.57 dB in Appendix A 
of [1] for the parameters of the Ku-band radar. Hughes also uses this 
value. 

The constant F.A.P. loss has been set at 1.7 dB. This is an 
approximation and is an area which deserves more study. 

The doppler mismatch is a loss of signal gain due to the fact 
that the doppler frequency is not at the center of one of the filters 
in the DFT. The loss varies from 0 dB to 3.9 dB, with an average of 
approximately 1.1 dB. 

For an N point DFT, the filter transfer function is given by [1]: 


F(Af) 



N sin (2TrAf -f) 


0 ) 


where Af = frequency difference from the center of the doppler filter 
to the doppler frequency of the target 
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Table 1. Summary of Target Detection Processing Losses 



Hughes 

(average) 

(dB) 

Axiomatix 
(worst case) 
(dB) 

Transmit Loss 

3,7 

3.7 

Beam Shape Loss 

2 (avg) 

3.2 (max) 

Pre-Sum Mismatch 

0,57 

0.57 

Constant FAP 

1.7 

1.7 

Doppler Mismatch 

1.1 (avg) 

3.9 (max) 

Range Gate Straddle 
Designated Mode 
Undesignated Mode 

1.76 

L str = f W 

1.76 

L s tr‘ f < R > 

Processor Loss 

1.25 

1.25 

T . Designated Mode 

'■ Ll Undesignated Mode 

12.08 

10.32 + L str 

16.08 

1A.S2 + L str 


Table 2. 

Doppler Mismatch Losses of DFT Filters 


N» 

1 

2 

4 

8 

16 

32 

Average Loss (dB) 

0 

-0.871 

-1.055 

-1.100 

- 1.111 

-1.113 

Maximum Loss (dB) 

0 

-3.010 

-3.698 

-3.867 

-3.908 

-3.919 
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Region of Pattern Overlap 



Edge of 3 dB Antenna Pattern 


Direction of Antenna Motion 


Target Location for Maximum 
Beam Shape Loss 


Figure 2. Antenna Beam Shape Loss Description 
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N = number of points in the DFT 
T = [PRF]~^ seconds. 

p 

The filters are placed (NTp) Hz apart, so that the maximum 
value of af is (2NTp)”\ at which value maximum signal strength moves 
to the adjacent filter in the DFT. 

The value of | F ( Af ) | at af=0 is 0 dB. The peak loss occurs at 
Af=(2NT )“ , which is given by 

r 


Maximum Loss 


|F(fif)| 2 



1 

Af 2NT p 

N sin ( ir/2N T 


( 2 ) 


The average loss is evaluated from 


Avg Loss 




sin (ux/2) 

N sin (irx/(2N)) 


2 

dx . 


(3) 


These doppler mismatch losses have been computed, and the results 
are shown in Table 2 for various values of N. The present radar has 
an N=16 point DFT, which results in 1.11 dB average loss and 3.90 dB 
maximum loss. These numbers are used for the doppler mismatch loss in 
Table 1. 

The range gate straddling loss [1, Appendix A] is 1.76 dB in the 
range designate mode . This is a factor of 1.5 (1.76 in dB), which arises 
since the range gates during acquisition are 1.5 times wider than the 
pulse width. In the undesignated mode , the range gate straddling loss 
is a function of range and is plotted in Figure 3. This variation in 
range is taken into account in evaluating detection performance. 

At a PRF of approximately 3 kbps, the unambiguous range is approx- 
imately 27 nmi. As the range gets close to multiples of 27 nmi, the SNR 
must get arbitrarily large for detection to take place, due to eclipsing 
of the received pulse by a new transmitted pulse. When the radar is 
operating in the cooperative mode with an active target , the unambiguous 
range increases to 300 nmi. The SNR clearly is a problem at ranges much 
above 10 nmi, and acquisition is not expected to be reliable. 
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Straddling Loss, 



Figure 3. Range Straddling Loss in Undesignated Mode of Detection 
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The last loss that has been accounted for is a processor algorithm 
loss Most of the algorithms are approximations to the ideal computation. 
Although the loss is small in each algorithm, the cumulative effect of 
all of these approximations is estimated to be approximately 1.25 dB. 

The totals of the losses in Figure 3 are used in evaluating 
radar detection performance in section 2.3. It is noted that the dif- 
ference between the average losses and worst case losses is 4.0 dB. In 
the next section, the differences in these losses and their implications 
are discussed. 

2.2 Radar Detection Performance Philosophy 

The losses enumerated in the last section demonstrate two approaches 
in establishing a radar detection performance philosophy. The differences 
center around what is to be taken into account when averages are performed. 
The "Hughes" approach is to average over all target positions and velo- 
cities when carrying out an average. The "Axiomatix" approach is some- 
what more pessimistic in that the average is earned out over thermal 
noise, target scintillation effects, and other radar parameters. This 
is done, however, for a target located at the worst possible velocity 
and worst possible angular position. 

We believe there is merit to both approaches. Therefore, the 
detection results to be subsequently presented are carried out under 
both performance philosophies. 

In the one case, the design philosophy is that 

Average P d = 0,99 (4) 

for a 1 square meter target. In the latter case, the design philosophy 
is 

Average > 0.99 (5) 

for a 1 square meter target for all target velocities and spatial 
positions. 

2.3 Radar Detection Performance 

Based on the losses discussed in the previous sections, the 
Ku-band radar detection performance is summarized. In Figures 4 and 5, 
the target detection design margins are plotted versus range for both 
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Figure 4. Target Detection Design Margin Versus Range for Worst 

Case Losses in Table 1 
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67 8 9 10 11 12 


Range (nmi) 

Figure 5. Target Detection Design Margins Versus Range for 

Average Losses in Table 1 
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the designated range and nondesignated range modes of operation. Except 

for the losses, all radar equation parameters are the same as used in 

[1, Appendix A]. In Figure 4, the worst case losses of Table 1 are 

used. Also, a transmitted power of 50 watts is used, although the 

present system configuration consists of a 60-watt TWT. This is 

approximately a 0.8 dB improvement. For the range designation mode, 

a design margin of 0 dB is obtained at approximately 9.3 nmi. For the 

nondesignated range mode, it is about 7.5 nmi. The design margin for 

the nondesignated range mode deteriorates as the range is reduced in 

the interval of 8 to 11 nmi because in this interval, the range gate 

-4 

straddling loss is increasing faster than the factor R in the radar 
equation. At 12 nmi, the design margin is -4.5 dB for the range desig- 
nated mode and -5.1 dB for the nondesignated range mode. 

In Figure 5, the design margin is shown versus range for the 
average losses in Table 1. In this case, a 0 dB design margin is 
obtained at approximately 11.7 nmi; at 12 nmi, the design margin 
is greater than -1 dB. 

The cumulative detection probability is plotted versus range for 
the worst case and average losses in Figures 6 and 7, respectively. 

Both M= 1 and M=2 scans are shown. It is seen that the two-scan case 
uniformly outperforms the single scan case for both the designated and 
nondesignated range modes. In these computations, a transmitted power 
of 60 watts was used. 

Inspection of Figures 6 and 7 reveals that, depending on the 
range, the cumulative detection probability for the designated range 
can be both smaller and larger than that for the nondesignated case. 

This is also the case for the design margins in Figures 4 and 5. At 
12 nmi, the cumulative detection probability is approximately 75% for 
two scans and approximately 50% for one scan This is the case for both 
the designated and nondesignated range modes. 

In order to determine what can be gained by allowing more total 
time to detect the target, the cumulative probability of detection was 
determined for various scan times and at various ranges. In Figure 8, 
the cumulative probability of detection versus total scan time is shown 
for various ranges. In these computations, two scans were assumed so 
that the scan rate varies with the total scan time. By allowing a slower 






Figure 6. Cumulative Detection 

Losses 
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Probability Versus Range for Worst Case 
in Table 1 





Cumulative Probability of Detection 
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Figure 7. Cumulative Probability of Detection Versus Range for 
Average Losses in Table 1 



Cumulative Probability of Detection 
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Figure 8. Cumulative Probability of Detection Versus Scan Time for 

Various Ranges 
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scan rate, it is seen in Figure 8 that there is dramatic improvement in 
the cumulative probability of detection. In particular, if the total 
scan time is increased to 2 minutes, then the cumulative detection 
probability at 12 nmi is 99% (the specification value). This is for 
the worst case losses and a designated range mode. 

Unfortunately, increasing the coherent detection time in this 
manner by lowering the scan rate has major hardware implications. It 
does provide the best improvement in performance, however. 

The scan rate is held fixed at two scans per minute in Figure 9, 
where the cumulative detection probability is plotted versus range. 

In this case, the improvement in performance is not as dramatic as when 
the scan rate is reduced. The gain is significant, however. At 11.3 nmi, 
for example, by increasing the total scan time from 1 minute to 2 minutes, 
the cumulative detection probability increases from 90% to 99%. 

2.4 Alternatives and Recommendations 

If the design philosophy of average losses is acceptable, then 
the design margin at 12 nmi is within a few tenths of Q dB with a 50-watt 
TWT and approximately 0 dB with a 60-watt TWT. 

If the design philosophy of worst case losses is preferred, the 
design margin is betv/een -4 dB and -5.2 dB. In this section, various 
alternatives are listed, along with recommendations. 

(a) Use a Slower Scan Rate . As described in the previous sec- 
tion, this provides the maximum improvement in performance for a given 
amount of increased time Since this alternative has major hardware 
implications, it is not recommended . 

(b) Increase Post-Detection Integration (PDI) . This can be 
accomplished by increasing the total number of RF frequencies used in 
frequency diversity This also has moderate hardware implications and 
is therefore also not recommended . 

(c) Increase the Number of Scans at the Same Scan Rate As shown 
above, there is significant improvement in performance by maintaining 
the same scan rate and allowing the total search time to increase from 

1 minute to 2 minutes. This increases the total number of scans from 

2 to 4. This keeps the present radar configuration entirely intact and 
thereby has virtually no hardware implications. For worst case losses. 



Cumulative Probability of Detection 
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Range (nmi) 


Figure 9. Cumulative Detection Probability Versus Range for Various Scans 
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a 0 dB design margin is obtained at a range of 11.3 nmi. He recommend 
this alternative , and this recommended design point is shown in Figure 9. 

(d) A final alternative which applies to the range designate 
mode only is to narrow somewhat the two range gates during search. How 
much they can be narrowed is a function of the anticipated range desig- 
nate accuracy. We do not have a recommendation on this alternative. 

It could be implemented in addition to the recommended third alternative 
above. If one has any doubt about the range designate accuracy, however, 
then this alternative does not appear satisfactory. 
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3.0 POINT TARGET TRACKING 


The tracking performance of a point target of the Ku-band radar 
is presented. The performance is described twofold. First, the effects 
of thermal noise and target amplitude scintillation are given in sec- 
tion 31. In section 3 2, the tracking capability is given in the pre- 
sence of relative accelerations along and perpendicular to the LOS. 
Emphasis is placed on the recovery time after an acceleration. 


3.1 Thermal Noise and Target Scintillation Effects on 
Tracking Accuracies for Passive Point Targets 


The purpose of this section is to present the numerical results 
associated with the analysis carried out in [1] for the Ku-band radar 
tracking accuracies due to thermal noise and target scintillation 
effects for a passive point target . The target effects which become 
a major concern at short ranges will be discussed in a subsequent sec- 
tion. No attempt will be made at this point to perform any additional 
analysis of the tracking loops. 

The tracking accuracies of the four variables being tracked 
(range, range rate, angle and angle rate) are presented as functions 
of range and compared to the existing specifications [2]. It is to be 
noted that there are two factors affecting the results: 

(1) Thermal noise and target scintillation which become the 
dominating factor at long ranges 

(2) Quantization effects in the Ku-band radar [1] which uses 
4 bits of quantization, resulting in a signal -to-noise. ratio at the 
output of the A/D converter equal to 


< SNR W 


1 

1 + (1 + (SNR) in ) (0.0129)_ 


($NR) in > 


( 6 ) 


where (SNR) in is the signal-to-noise ratio at the input of the A/D con- 
verter. The signal at the output of the A/D converter is presumed with 
a resulting gain in signal-to-noise ratio equal to 16 (12.04 dB). Thus, 
the resulting SNR used in the calculations is 

(SNR) dB - «SNR) out ) dB + 12.04 dB. 


(7) 
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It is easy to see from (6) that, at high input signal -to-noise ratio 
(close range), the resulting SNR saturates at 

SNR = 18 9 + 12.04 = 30.94 dB , 

which makes the quantization effect the major contributing factor to 
the tracking accuracies. The combined effect of both factors is taken 
into account in all subsequent calculations. 

There are many parameters in the various tracking loops that are 
varied with the designated range. The parameters used in these calcu- 
lations are the latest available from Hughes Aircraft Company and are 
summarized in Table 3. 


3.1.1 Angle and Angle Rate Tracking Accuracies 

The angle and angle rate tracking accuracies are presented together 

because the angle rate estimation is carried out in the angle tracking 

loop. The ratio of angle rate root-mean-squared (RMS) error a • to the 

angle RMS error a. is given as [1]* 

0 



Kt 2 + 1 



(8) 


where K is the equivalent angle loop gam and t is its time constant. 
The variations of K and t with range are shown in Table 3. The angle 
RMS error is given by 


(Kt £ + 1)T 


4t 


tA: °x9< n ’ SNR > 


9 m 


(9) 


where T g = time per sample of the equivalent discrete time noise (x) 
whose variance is a a 

X O 

e g = 1.6° = 3 dB bandwidth of the Ku-band antenna 

v 1 

N " number of independent samples used to obtain a „ (N = 1 0 

XU 

for angle tracking) 

p. and cr . are functions of SNR and are discussed in detail in 

□ XU 

[1] and in the simulation section of this report. 


Figures 10 and 11 illustrate the variations of 3c. and 3a- with 

U u 

range. It is obvious that both the angle and the angle rate accuracies 



Table 3. Summary of Tracking Loop Parameters (Passive Target) 


Parameter 


>9.5 

0.099 


Range (nmi) 


9. 5-3.8 3.8-1 9 1.9-0.95 0.95-0.42 


: 1 

i 

(msec) j 

5.36 

A (usee) 

33.2 

a | 

0.0566 

- — ■ — j 

8 

0.000884 

Bp (Hz) 

186 

m 

2 

T $ (msec) 

107.2 

L e,0 

16.90 

1 

1 . 
l R,R 

9.49 

K 

0.0288 

T (sec) 

12 

PRF (Hz) 

2987 


0.116 


2.29 


16.6 


0.1132 


45.9 


16.83 


9 42 


0.0288 


12 


0 058 


2.29 


8.3 


0.2263 


0.00354 0.0283 

437 




45.9 


16.83 


9.42 


0.2221 


4.25 


697 


0.029 


2.92 


4.15 


0.2263 


0.0283 


437 


4 


45.9 


16.83 


9.42 


0.5685 


0.0144 


2.29 


2.07 


0.4526 


0.1132 


437 


45.9 


16.83 


9.42 


0.5685 


Remarks 


0.00085 1 Duty factor 




0.122 Pulse width 


0.4803 


0.1202 


Doppler filter 
bandwidth 


Number of 
averaged samples 
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9.42 


0.5685 
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6970 




Angle and angle 
rate losses 
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rate losses 


Angle loop gain 


Angle loop 
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Figure 10. Angle Tracking Error of Passive Point Target Due to Thermal 
Noise and Target Scintillation Versus Range 



Angle Rate, 3a* (r.irad/sec) 
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meet the design specifications. The design margin in the angle rate 
tracking, however, is much less than that for the angle. This is due 
to the fact that the loop filter in the former has been widened to pro- 
vide faster recovery from acceleration transients. It is believed that 
increasing the gain K slightly would improve the transient performance 
of the loop while maintaining the random error within or very close to 
the specifications. 

3.1.2 Range Tracking Accuracies 

The RMS of the range error o ^ (timing error) is computed using 
the results of [2], namely, 

1 

°R ■ TV + sSr) (B N T s ) 1/2 o xR (N,SNR), (10) 

where A = pulse width of the transmitted signal (see Table 3) 

c = velocity of light in free space 
= bandwidth of the range tracking loop which includes the 
a - 3 tracker as a loop filter 

B . aU+VSHRp'y + B („) 

N 4aT s ' ' 

a xR = ec l u ' ,va l ent noise of the system due““t‘o thermal 

noise 

N = number of independent samples used to obtain a x ^ (N = 20 
for range tracking. 

Figure 12 illustrates the variation of the RMS range error versus 
range with the noise bandwidth as a parameter. The gump at 9.5 nmi is 
due mainly to the change in PRF at that range. It is observed that, 
since the values of a and 3 used at longer ranges correspond to £5^= 1 Hz, 
the range RMS error due to thermal noise meets the specifications up 
to 9.5 nmi. The specifications are illustrated in Figure 13 for 
convenience 

Figure 14 illustrates the variations of the range error at closer 
ranges (R<6 nmi) with the latest values of c and 3 as parameters. The 
figure also shows the latest performance curve presented by Hughes for 
comparison purposes. There are four points to observe. 
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Figure 12 RMS Range Tracking Error of Point Passive Targets Versus Range 
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Figure 13. Specification for Three-Sigma Random Range Error Versus LOS Range 
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Figure 14. RMS Range Tracking Error of Passive Point Targets at Close Range 
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(1) The range tracking performance of a passive point target 
due to thermal noise and target scintillation is well below the speci- 
fications, especially at close range. 

(2) The flat portions of the curve at ranges R<1 nmi are due 
to quantization limiting effect. 

(3) The difference between the Axiomatix and Hughes performance 
curves is believed to be due to the automatic gain control (AGC) which 
has not been taken into account in the analysis. 

(4) Although the accuracies at close range drop below 10 feet, 
it is believed that 10 feet is the best achievable accuracy because of 
the sensitivity of the meters. 


3.1.3 Range Rate Tracking Accuracy 

The range rate RMS error due to thermal noise and target scintil- 
lation has been found in [2] to be 




(N.SNR) , 


( 12 ) 


where Bp = doppler filter single-sided noise bandwidth 

f - RF carrier frequency. The middle RF frequency 

f = 13.885 GHz in the frequency diversity system is 
used in the computations. 

m - number of samples being averaged to give a single R 
reading. m= 2 for R > 9.5 nmi and m - 4 for R< 9.5 nmi. 
a ^ = RMS of the equivalent noise of the system due to thermal 
noise. 

N = number of independent samples used to obtain (N=40 
for range rate tracking). 

The evaluation of for various ranges is shown in Figures 15 
and 16. Figure 15 illustrates the range rate errors at long ranges and 
shows the effect of varying the number of averaged samples (m) at 9.5 nmi. 
As is obvious from the figure, the range rate accuracies meet the speci- 
fications for all ranges below the specified 10 nmi. The same thing is 
true for Figure 16, which illustrates the variations of at short 
ranges. 
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Figure 16. RMS Range Rate Tracking Error Due to Thermal Noise and Target Scintillation 
of Passive Point Targets Versus Range (Short Ranges) 
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In order to avoid confusion, it should be pointed out that the 
numerical evaluations of the range rate accuracies in [1] were erroneous 
and the numbers presented in Figures 15 and 16 represent the correct 
calculations. 

3 2 Effects of Point Target Accelerations 

3.2.1 Introduction 

This section investigates the performance of the Ku-band radar 
system in the tracking mode when a constant acceleration takes place. 

Two extreme cases can be distinguished. The first will be referred to 
as radial acceleration , which is achieved by firing retro rockets m 
the direction of the line-of-sight (LOS). This causes the range to vary 
as a quadratic function of time rather than a linear function of time, 
which occurs in the case of a constant approaching velocity (range rate). 
The second case will be referred to as tangential acceleration , which 
is an acceleration in a direction perpendicular to the LOS, causing a 
quadratic variation in the target angle being tracked. Since both the 
azimuth and the elevation angles are tracked similarly [1], this section 
will not distinguish between the two angles. The radial and tangential 
acceleration cases are illustrated in Figure 17. 

The analysis of the tracking loops in the presence of accelera- 
tion is performed in this section for a deterministic input. This is 
a realistic assumption for the range of interest, namely, close range 
when the signal-to-noise ratio is high. It is believed that the tran- 
sient behavior of the tracking loops with noisy inputs can only be 
determined via simulation of these loops, which is discussed in Section 
4.0 of this report. 

3.2.2 Radial Accelerations 

A constant radial acceleration, as shown in Figure 17, causes the 
range to change as a quadratic function of time without affecting the 
angle measurements. The range tracking loop is discussed in [1] and is 
shown in Figure 18 for the deterministic case of interest ($NR-*-“> implies 
p=l and no noise). The parameters in the figure are: 

T = actual time delay between transmission and reception of a 
given pulse. For a passive target at range R, T=2R/c, 
where c is the velocity of light in free space. 


QMGWALP*®® 
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Radial 
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Figure 17. Radial and Tangential Accelerations 
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Figure 18. Deterministic Range Tracking Loop 
T = loop estimate of the actual delay T. 

A 

Ay = I - T = the timing error. 

T = computation time per sample = 107.2 msec when R>9.5 nmi 
s = 45 9 msec when R<9.5 nmi. 

a, 8 = a-8 tracker parameters which vary with the designated range 
as shown in Table 4. 


Table 4. a-8 Tracker Parameters 





Range (nmi) 



<0.42 

<0.95 

<3 8 

<9.5 

<13.2 

a 

0.4809 

0.4526 

0.2263 

0.1132 

0.0566 

8 

0.1202 

0.1132 

0.0283 

0.00354 

0.000884 


Simple derivation results in 

s 2 ORIGINAL PAGE IS 

AT(s) = — g -ryr rgyr T(s) , OF POOR QUALITY 

s +8 lV s+ (? r ) 

s 

which, when converted to range measurements, yields 


s^ + a‘s+ y— 


AR(s) 


R(s) , 


(13) 
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where R(s) = Laplace transform of the actual range R(t) 
AR(s) = Laplace transform of the range error AR(t) 

a' and 3 ' are defined as a* = ^ and 3 ' = ^ . 


Since the actual acceleration (A) is a rocket burn for a period 
of time T, it can be represented as 

a{t) = A [u ( t) - u(t-T)] , (14) 


where u(t) is the unit step function. Taking the Laplace transform of 
(14) and using the fact that the acceleration is the second derivative 
of range, the range transform can be written as [3]: 

R(s) = 4-(l-e" sT ) (15) 
s d 


which, when substituted in (14), results in: 


AR(s) 


A(1 - e~ sT ) 


S 2 + a* 


s+ V 


(16) 


The denominator has a positive discriminant for all the values of a and 
3 listed in Table 4 and hence can be factored as sts + a-jHs + ag), where 

a-j - y ^CC - a - 3/2 ^ 

a 2 ■ f (“+ 7 « 2 - e/2 ) 

S x 

Taking the inverse Laplace transform of (16) yields 


AR(t) = A 


1 “ 2 ' 
+ 




-a 2 t\ 


a l a 2 


a l a 2 {a l - a 2> 


u(t) 


1 +^1 


■a, (t-T) -a ? (t-T)\ 

* - a a (- 


- a-j e 


a l a 2 


a l a 2 (a,-a 2 ) 


lu(t-T) 


(17) 
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Figure 19 illustrates the values of AR as a function of time at 

2 

a range of 10 nmi with an acceleration A=0 Ig (g=32 ft/sec ) and a 
burn time T- 10 sec, while Figure 20 shows the variations of AR for 
A = O.lg and T = 10 sec at a range of 2 nmi. It can be easily seen that 
the error due to O.lg acceleration is well below the la specifications 
for the range error (lc=40 ft and 100 ft at 2 nmi and 10 nmi, respec- 
tively). Since it is believed that the value of the acceleration (A) 
will actually be of the same order of magnitude as O.lg or less, no 
problems are anticipated in tracking radial accelerations during target 
approach. 


3.2.3 Tangential Acceleration 

While it was shown that the radial acceleration affects the range 
measurements, the tangential acceleration has a direct effect on the 
angle and angle rate measurements. The deterministic angle and angle 
rate tracking loop [1] is shown in Figure 21. In the figure. 


e T = actual target angle 

e A = antenna angle 

A e = 9j- 0^ = angular error 

A 

Qj = estimate of the angle rate 
K = effective loop gain 


x = time constant 


K and t are given in [4] as a function of range and are listed 
in Table 5. 


Table 5. Angle Tracking Loop Parameters 


Range (nmi) 

<1.9 1.9-3 8 >3 8 

K 0 5685 0.2221 0.0288 

t 2.7 4 25 12 






AR (ft) 



Time (sec) 


Figure 20. Range Error Due to Constant Acceleration Lasting T Seconds 

(R = 2 nmi ) 
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Figure 21. Deterministic Angle Tracking Loop 


Using Figure 21, it can be written that 

A (s) = e T (s) 

e s 2 + K(st + 1) T 


L ( s) = Ks e < s ). 

T s 2 +K(st + 1) T 


(18a) 

(18b) 


However, assuming that the retro rockets will cause a constant tangen- 
tial acceleration A for T sec, the target angle transform can be written 
as 


e T (s) 



±-0-e- sT ) 


(19) 


where the range R is assumed much larger than the tangential displace- 
ment that takes place in T sec. 

Combining (18) and (19) and taking the inverse Laplace transform 
yields the following results for the angle error and the angle rate 
error. 
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A 

R 


Ae(t) = 


—3-t t “3o"t 

1 , a 2 e 1 - a 1 6 


a l a 2 a 1 a 2 (a^ - a 2 ) 


; t < T 


( 20 ) 


Ra ] a 2 (a-, - a 2 ) 


a-iT -apt a 2 T , 
a 2 e (l-s ) - e (1-6 ) 


; t >.T 


KA 

R~a-j a 2 (a-j - a' 2 )' 


‘"Cli t 

t(a r a 2 ) + 57 (, - e > 

t < T 


Ae{t) = 


. A t . 

R z * 


(21) 


AK 

Ra iV a i- a 2 ) 


a-, -a 9 t a 0 T a 0 -a,T a,T 
T < a l- a 2> + 4 e 2 0-eZ)-^e 1 (1-e 1 ) 


_A_ . 

R * 


t>T 


where 


>i - H kt -J k2t2 - 4i< ) 

> 2 ■ K Kt+ 7 k2t2 - 4k )- 


As is the case in the range tracking loop, the system is slightly over- 
damped for all values of K and t given in Table 5. (The discriminant 
of the denominator of (18) is positive.) 

Figure 22 illustrates the behavior of the angle error for a burn 
time T=10 sec and an acceleration A=0.1g for various ranges, while 
Figure 23 illustrates the angle error at a range R= 1000 ft for a 
number of burn times. Comparing the obtained results with the angle 
random error specifications, it can easily be seen that the error due 
to acceleration is always less than 3cr and is less than ler for all the 
cases considered except for ranges less than 1000 ft and burn times 
larger than 2 sec. In the latter two cases, the recovery time, which 
can be defined as the time required after the burn to reside below 
the specifications (ler), is equal to 4 sec. Since the actual magnitudes 



Ae(t) (mrad) 
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1 R = 10 nmi ; f = 0.027 Hz; 2x - 24 sec; = 0.107 Hz 

2 R = 2 nmi ; f n = 0.075 Hz; 2x = 8.5 sec; = 0.29 Hz 

3 R = 1000 ft; f = 0.12 Hz; 2t = 5.4 sec; b‘ n - 0.48 Hz 



Figure 22. Angle Error Due to Tangential Acceleration With Range as 

a Parameter 



ke(t) (rirad) 
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of acceleration are believed to be in the range of O.lg or less, the 
angle errors due to tangential accelerations are not believed to con- 
stitute a major problem for angle tracking. 

Similar to the angle error variations, the angle rate error vari- 
ations are illustrated in Figures 24 and 25 for the same parameters as 
used in Figures 22 and 23, respectively. It is evident that the angle 
rate errors due to acceleration exceed the 3a angle rate specifications 
for all cases studied. For the close range of 1000 ft, the recovery 
times to reside within la of the angle rate specifications are 8, 8.75, 
9.25 and 9.30 sec for burn times of 1 , 2, 5 and 10 sec, respectively. 

All of these recovery times exceed the desired 2 sec recovery time . 

In order to reduce the angle rate acceleration errors to the desired 
levels, the equivalent noise bandwidth of the angle tracking loop has 
to be widened considerably, which would result in an appreciable 
increase in the angle and angle rate errors due to thermal noise. It 
is believed that, short of major design changes, any attempt to reduce 
the angle rate recovery times to within the desired value of 2 sec is 
unrealistic. It is therefore recommended that no design changes be 
considered due to the fact that angle rate estimation does not meet 
specification values within 2 sec after a burn. 

Responses from other values of accelerations for any of the 
cases in this section can be obtained simply by scaling the desired 
value with respect to that shown on the appropriate figure. This is 
due to the fact that, at high values of SNR and the small accelerations 
considered, the systems are nearly linear. The same scaling cannot be 
carried out in the time dimension, however. 

3.24 Angular Acceleration Testing 

One of the radar testing procedures proposed by NASA is to fly 
a target (airplane) at an altitude R in a circular path of radius r 
and an angular velocity of to radians/sec, as shown in Figure 26. 

The angular acceleration being tracked by the radar is then* 

2 

e'-j-(t) = — cos tot . (22) 

ORIGINAL PAGE IS 
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1 R = 10 nmi ; f = 0.027 Hz; 2 t * 24 sec; = 0.107 Hz 

2 R = 2 nmi; f = 0.075 Hz; 2 t = 8.5 sec; = 0.29 Hz 

3 R = 1000 ft; f = 0.12 Hz; 2x = 5.4 sec; « 0.48 .Hz 



Figure 24. Angle Rate Error Due to Acceleration With Range as a Parameter 



Ac(t) (mrad/sec) 
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Figure 26. Radar Angle Tracking Testing 


The Laplace transform of the angle being tracked can be simply shown 
to be [3]: 


9j(s) = 


rw 


Rs (s 2 + to 2 ) 


Substituting (23} in (18) results in 

Ms) = 


„ 2 

r to s 


R(s 2 +co 2 )(s 2 +K(ts + 1)) 


and 


e T (s) 


r (/ K 


R(s 2 + to 2 ) (s 2 + K (ts + 1)) 


(23) 


(24a) 


(24b) 


Taking the inverse Laplace transform of (24a), and after considerable 
algebraic manipulation, the angle and angle rate errors are found to 
be 


■ nr 


ye 


-yt 


fie 


-fit 


(y-6) (y 2 + (i) 2 } (fi~y) (fi 2 + W 2 ) / (y 2 + U 2 ) (fi 2 + 03 2 ) 


COS (oit+lp) 
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The first two terms represent the transient behavior; hence, they can 
be dropped m the steady-state response calculations, which results in 

A (t) = " (y 2 + to 2 )’' 1/2 C<5 2 + to 2 )“ 1/2 cos (ut+*) , (25) 

e ss * 

where 

+ ■ t + ta,rl (£) - tan ' 1 (t) 



As expected, the steady-state response is a sinusoidal function of time, 
as is the input with a different amplitude and phase shift. 

The steady state of the angle rate response can be derived simi- 
larly, yielding 


A* (t) 
e ss 


A / (y6-<h Z - K) 2 + tt 2 (S + Yr f + . 

“ ' (TTTuTTTi ’ 


(26) 


where 


<i> = tan‘ 


>(s + y) 


[j 6 - to - Kj 


+ tan_1 (£) - tan_1 (t) 


Equations (25) and (26) provide the required expressions for the theor- 
etical evaluation of the angle and angle rate performance in the proposed 
testing and thus establish the basis for comparing the experimental 
results with the theoretical analysis. 
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4.0 SIMULATION AT SHORT RANGE 

The results of the previous sections, target effects and other 
effects are combined in this section to develop an overall model of 
the Ku-band radar at short ranges (less than a few thousand feet). The 
extension of the model to longer ranges does not involve extensive modi- 
fication. This model is developed so that it can be programmed directly 
into a digital computer simulation. 

There are various simulations at NASA JSC of the Orbiter and, in 
particular, of the proximity and docking operations. Some have a man in 
the loop and others do not. Some operate real-time and others do not. 
The simulation model presented here is developed so that it is essen- 
tially modular. That is, with minor interfacing, the simulation model 
can be used in any of the simulations presently operating at NASA JSC. 

The simulation model will interface with the existing simulations 
as shown in Figure 27. The existing Orbiter simulations continuously 
provide a "trajectory information vector" to the radar simulation. 

Among the components of the trajectory information vector, T(t), are 


T(t) = 


R<t) 

R(t) 

R(t) 

0 AZ^ 

e AZ (t) 

0 AZ (t) 

0 EL^ 

e EL (t) 

e EL (t) 

L R (t) 

L 0 (t) 


Range along LOS 

Velocity along LOS 

Acceleration along LOS 

Azimuth angle 

Azimuth angular rate 

Azimuth angular acceleration (27) 

Elevation angle 

Elevation angular rate 

Elevation angular acceleration 

Effective target length along LOS 

Effective target length perpendicular 
to LOS 


The components of T(t) in (27) are not intended to be complete. 

The radar simulation assumes these components are available when necessary. 
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Figure 27. General Block Diagram of Orbital Simulation and the Simulation of 
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Other necessary inputs to the radar simulation are internal system 
constants, e.g., PRF, pulse width, and a set of mutually independent 
random noise generators. Both uniform and Gaussian random noise gener- 
ators will be needed. 

The radar simulation will output an error vector, each component 
of which will represent a tracking or an estimation error. Each of 
these errors will be the result of several error contributors. The 
effect of each contribution will also be available. The error vector 
is sent back to the existing simulator when some or all of the error 
components are added to the components of T(t). The measured or observed 
variables are then given by 

0 ( t) = T(t) + E(t) . (28) 


In the succeeding parts of this section, the necessary components 
of the radar simulation are detailed. Finally, the overall description 
is examined again at the conclusion of this section. 

4.1 Simulation of Thermal Noise and Target Scintillation Effects 


In order to simulate the errors due to thermal noise and target 
scintillation on a digital computer, discrete models fo r th e tracking 
loop operations must be developed. In the Ku-band radar, there are 
two primary closed tracking loops to track range and angle. The range 
rate is tracked via an open loop , and the angle rate is estimated in the 
angle tracking loop [1]. 


4.1.1 Tracking Error Simulation Due to Thermal Noise and Target 
Scintillation in Range, Angle, and Angle Rate Tracking 


The range tracking loop is shown in Figure 28, while the angle 
tracking loop is shown in Figure 29. All the parameters are defined 
in [1, Appendices E and F] and will not be repeated here. 

In the range tracking loop, the closed loop transfer function 
between the noise input n and the estimated time of signal arrival 
(range) T is given by [1, Appendix F]. 


H r (s) 


a' (S+B7(«' T s )) 
s 2 +a's + B'/T s 


(29) 



— 1 > 
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where a' and 3 ' are related to the a- 3 tracker parameters via the equa- 
tions: 


A 8p R e 


T is the time per sample of the Gaussian random noise n v by the D/A 
converter whose mean is zero and whose variance is equal to [1]. 
This variance is a function of the signal-to-noise ratio and the number 
of independent samples (N) being post- integrated. For range tracking, 
N= 20. 

Similarly, for the angle tracking loop, the closed loop transfer 
function between the input noise n Q and the angle estimate is given 
by [1 , Appendix D] 


H e {s) 


K(ts+1) 
s^ + K(tS + 1 ) 


(30) 


whereas the associated transfer function for the angle rate is [1, 
Appendix D]: 


H j(s) 


Ks 

s^ + K(ts + 1 ) 


(31) 


where K is the equivalent loop gam and t is a time constant defined 
in [1], 

The noise n Q is Gaussian random noise with zero mean and variance 
equal to 



2 

The variance ct x0 , as is the case in the range tracking loop, is a func- 
tion of the signal-to-noise ratio and the number of independent samples 
(N) being post-integrated. For angle tracking, N=10. 

It is evident from (29) through (31) that the three transfer 
functions can be parameterized using the general formula 
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a 4 ( s +aJ 

’ (s+a 1 )(s + a 2 ) * 


(33) 


This is possible since for all the parameters under consideration the 
denominators of the transfer functions have positive discriminants. 

The latest values of a, 3, t, and K are summarized in Table 6, while 
the parameters a^, a^, a 3 , and a^ are defined for the three cases in 
Table 7. It is noted that a-j , a 2 » a 3> and a^ are all positive numbers 
and that a 1 > a^. 

4.1.2 Stochastic Process Modeling 

The continuous stochastic process modeling of the continuous 
tracking loops can now be obtained as shown in Figure 30. The equations 
describing the processes are 

x-|(t) = -a 2 x 1 (t) - (a 1 > a 3 ) x 2 (t) + a 4 u(t) 

x 2 (t) = -a ] x 2 (t) + a^u(t) , 

which can be written in a matrix form as 

X{t) = A X(t) + B u(t) , (34) 

where 


A 

X-| ( t) 

, A = 

“ a 2 ” (a r a 3 } 

• i i 

a 4 


“ l 

4^ 

CM 

X 

J 


0 -a.| 


_ a 4_ 


u(t) 


a 4 

x (t) 

s + a 3 

X-, (t) 


s + a-. 

i 

tr 

s + a 2 



Figure 30. Continuous Stochastic Model 
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Table 6. Tracking Loop Parameters 



Range (nmi) 

Parameter 

<0.42 

0.42-0 95 

0.95-1.9 

1.9-3. 8 

3. 8-9.5 

>9 5 

a 

0.4809 

0.4526 

0.2263 

0.2263 

0.1132 

0.0566 

e 

0.1202 

0.1132 

0.0283 

0.0283 

0.00354 

0.000884 

K (rad) 2 

0.5685 

0.5685 

0.5685 

0.2221 

0.0288 

0.0288 

(sec) 

2.7 

2.7 

2.7 

4.25 

12 

12 

T s (msec) 

45.9 

45.9 

45.9 

45 9 

45.9 

107.2 


Table 7. General Transfer Function Parameters 


a d (s+aj 
” (s + a 1 ){s + a 2 ) 


Loop 

a l 

a 2 

a 3 

a 4 

Range 

o' +J a ,2 -407T s 

BW 

8' 

a' 

2 

2 

ct'T 

s 

Angle 


| hh 

2_ 

Kt 

Kt+Jk 2 t 2 - 4 

SfefS 

2 

2 

X 

Angle Rate 


mmsm 

0 

K 

2 

2 
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In (34), the parameter being tracked is x,(t) [R, e or e]. Also, 

2 

u(t) is white Gaussian noise with zero mean and a variance a equal to 

O 2 u 

o D for the range tracking loop and equal to o„ n for the angle (angle 
xk ne 

rate) tracking loop. 

The discrete time model of (34) can be written as 


X(k+1) = CX(k) + Du(k). (35) 

In Appendix A of this report and in [5], matrices C and £ are given by 

C = e- T (36) 

F = D[Eu(k)u'(k)j D' = of fe^BB' e^- x dx , (37) 

u Jq 


where T is the computer update time. 

A T 

There are numerous methods for finding e— [5-7]. One simple 
approach is 

e- T = ^ _1 {[sl- A]" 1 } , (38) 

where is the inverse Laplace transform. Performing (38) yields 


C = 





£ 




where 




(39) 


(40) 


To calculate £ and the covariance matrix of u.(k), namely, [E(u_(k) u_'(k))], 
the following integration has to be performed 


F 


( a 4 ff u) 


•0 


-a ? x 
e L 



/ -a-,x\ 

e(e c - e ) 

i i 

-a ? x 

e 

0 

-a^x 

e 

i i 

( -a 2 x ny 
e( e - e 

) T* 1 * 


dx 


which reduces to 
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Since A is in Jordan canonical form [7-8] and thus cannot be diagonalized, 
or since x-j(t) and x 2 (t) are not independent, £(k) is a two-dimensional 
discrete process which implies that two Gaussian random generators are 
needed to realize this discrete process. 

Assuming that the variance of j^(k) is given by 

” 2 
a l 

E[u(k)u'(k)] = 

0 



which implies that the two components of u(k) are independent and 
expressing the matrix £ as 



then 

F = D E(_u(k) u!(k)) D_J 


j 2 2.2 2 

d n 0 i +d i2 a 2 


d ll d 21 a l 2+d 22 d l 2 a Z 


d ll d 21 0 l 2+d 22 d 12 a 2 2 
d 21 0 ] 2+ci 2 2 2 0 2 2 


( 42 ) 
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Equating (41) and (42) results in three equations to be solved 
for the entries of £ and E (u.(k) u/ (k)) . However, since there are six 
unknowns involved, three have to be chosen arbitrarily without violating 
the physical nature of the processes. Setting a-j =1, d-j^ = l, d22 = 0» 
and solving for the remaining entries results in: 


1 -e 


-2& } T 


j - 2 1 

21 7 -(a.+aJTv 7 ^2a7T 

(1+e)(l-e 1 2 )\ e(1_,e 1 


d-j ^ 3^ 


2a 


1 


(43) 


a o = 


(l+s)(i-e 


•(a^jT 


a 1 a^ 


-2a, T -| 

I - eOr-e...... ' ) 


2a, 


1 -e 


-2a-jT 




2a, 


(44) 


and 

4 


-2a 9 T 9 

(1-e 2 )(1 +e) 2 

2a„ 


-(a,+a 9 )T 9 -2a,T“ 

2e(l+e)(l-e 1 2 ) . e 2 (l-e 1 ) 

*r ' 0“ ^ 


a l +a 2 


2a, 


, ,2 2 

( a 4 0 u "°2 

(45) 


A diagram of the discrete time model is shown in Figure 31, which results 
from (35) with d^ = 1, d 22 = 0, and = 1 . 

It is emphasized that, in order to simulate the thermal noise 
and amplitude scintillation effects for the range and angle tracking 
loops, two independent Gaussian noise generators are needed 

4.1.3 Software Implementation 

The actual implementation requires a software equivalent of 

Figure 31, where {a , I s5 1,4} can be calculated for a given range using 

1 2 

Tables 6 and 7. The parameters d 2 i , d-j 2 » and o are calculated using 

(43) through (45). In order to compute dg-j, d^, and o 2 , however, a ^ 

(the variance of the input noise in the continuous model) has to be cal- 

2 

culated. As stated previously, a is given by 
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' 6 B \ 2 



for range tracking 


(46) 


for angle and angle rate tracking 


2 2 

where a x0 and <? xR have been calculated in [1] and are shown in Figure 32 

as a function of signal-to-noise ratio (SNR) which is obtained from the 

o 

radar equation- Figure 32 also shows cr^, which is the variance of the 
input noise m range rate tracking, which will be discussed next. SNR 
is defined as the signal -to-noise ratio at the output of the doppler 
filters, SNR - P QV g r /{NgBp). For the Ku-band radar antenna, 0g=1.6° 
and k m =l. The only two additional parameters needed are the normalized 
correlation coefficients which are given as 


P R 


1 



(47a) 


for range, and 


1 + 


SNR 


(47b) 


for angle and angle rate. 

4.1.4 Tracking Errors Simulation Due to Thermal Noise and 

Target Scintillation in Range Rate Tracking 

The simulation of the thermal noise and target scintillation 
errors m tracking range rate is simpler than the previous cases, due 
to the fact that range rate tracking is done via an open loop, where 
m successive readings are averaged to give the desired range rate esti- 
mate [1, Appendix G]. This is accomplished via a sliding window with 
m-2 for R>9.5 nmi and m=4 for R<9.5 nmi. Contrary to range tracking, 
the range rate estimate is not used to continuously update the locations 
of the 32 existing doppler filters but rather to select the doppler 
filter closest to the range rate reading [1]. To carry out the range 
rate simulation process, it is first necessary to compute the signal - 
to-noise ratio via the radar equation and then to compute the normal- 
ized correlation coefficient which is given by [1, Appendix G]: 


original page IS 
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Figure 32. a x Versus SNR for Range, Angle and Range Rate Tracking Loops [1] 
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P R = 


(48) 


1 + 


•FT 


SNR 


The variance cr^ which is shown in Figure 32 is then computed, 
and the range rate mean square error per sample (<fy s ) is found using 
the following relation derived in [1, Appendix G]: 


c B r 


T R/s 


J 6 P R f Cy 


xR 


(49) 


where Bp is the doppler filter single-sided noise bandwidth, c is the 
velocity of light in free space, and f is the Ku-band RF frequency used. 
A white Gaussian noise generator is used with zero mean and a variance 
given by (49) to generate the required samples, and a sliding window 
averaging of the last m samples is used to give the range rate error, 
as illustrated in Figure 33. 

The simulated errors in range, range rate, angle and angle rate 
tracking due to thermal noise and target scintillation are combined with 
the other types of error in each case of interest to produce the overall 
simulated tracking error. 

4.2 Radar Simulation of Relative Accelerations 

The responses of the various Ku-band radar tracking loops to 
accelerations orthogonal to and along the line-of-sight (LOS) have been 
analytically developed in section 3.2 of this report. The results of 
numerical computation are also presented in section 3.2. In this sec- 
tion, we present digital simulation models which allow for arbitrary 
accelerations. When simulated real time, this allows for evaluating 
acceleration effects when a human operator is in the loop. 

4.2.1 Closed Loop Transfer Functions 

In the absence of noise, the range tracking loop is shown in 
Figure 34. Using this figure, the closed loop transfer function between 
the range R and the range error is (in the complex frequency domain): 



AR(s) 


R(s) . 


(50) 
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In terms of acceleration along the LOS, (50) becomes 


AR(s) 



(51) 


We therefore define the transfer function between acceleration along 
the LOS and range error as 


H r (s) 



(52) 


For the angle tracking loop, in the absence of noise, the block 
diagram reduces to that in Figure 35. The relationship between the 
target angle e^. and the angle error A0 is 

e T (s) , (53) 


A6(s) 


IS + K(tS+ 1 M 


which, in terms of the angular acceleration, becomes 


Ae(s) 


1 

s 2 + K(ts+ 1 )_ 


e' T (s) . 


(54) 


The closed-loop transfer function between angular error and angular 
acceleration (J_ LOS) is defined as 


H e < s > 


1 

s 2 + K(ts + 1 ) 


(55) 


Finally, also from Figure 35 for angle rate, the relationship between 
the angle rate estimate and the angle rate 0^. is 

Ks „ / _ _ K 


0-|-(s) p 

s + K(st + 1) 


®t(s) “ o ®-r(s) 

1 s^ + K(st + 1) T 


(56) 


. a • • 

In terms of the error of the angular rate estimate, A0 = - 0 T , and 

angular acceleration, the relationship is 
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from which we define the transfer function relation angular acceleration 
and angle rate estimate error as 


«e( s ) 


S + Kt 

s 2 + K(st+ 1) 


(58) 


For acceleration, the three transfer functions of primary interest 
are given in (52), (55), and (58). We note that H R (s) and H (s) are 
of the form 


H l (s) = Ts+'^Hs + aj) < 59a > 

and H*(s) is of the form 

^2^ “ (s + a 1 )(s + a 2 ) * (59b) 

The discriminants of the transfer functions are shown in Table 8, 
where it is observed that a ] and a 2 are real and positive in all cases. 
Digital simulation algorithms will therefore be developed for H-|(s) and 
H 2 (s) in (59) under the assumption that a > 0 Also note that a, >a 0 
in all cases. 


Table 8. Discriminants for Acceleration Transfer Functions 



Range 

Angle 

Angle Rate 

a l 

4 1 
V 


. — \ 


1 

2 

Kt + ^ 





U 2 - 6/2^ 

^K 2 t 2 - 4K | 

1 

2 


^K 2 t 2 - 4K 

a 2 

J-l 

V 

\ 

^a 2 - 8/2 j 


1 

2 




1 

2 



/|< 2 t 2 - 4K 

k 2 r 2 - 4kJ 

a 3 


— 


Kt 
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4.2.2 Digital Simulation Algorithms for Acceleration 

The inverse Laplace transforms of H^s) and H 2 (s) are 

= (a ] - a 2 ) ^ exp ^" a 2 t ^ " exp (60a) 

h 2^ = (irprii^)^ a 3 “ a 2^ ex P (— a 2^) — ( a 3 “ a i ) ex P (- a ]t)] . (60b) 
The z- transform of h^t), i = l,2, is given by [9-10]: 

03 

H,U) = I h.(nT)z“ n , i = 1 ,2 . (61) 

n=0 


For i-l, after some algebraic simplification, the z-transform 
can be written in closed form as 


H-|(z) = 


1 


ia r a 2 


2/ ~( a i +a 2^ T 


z V* T -e^ T ) 


z ~(e 


) . z^pV^) + lj 


(62) 


For i = 2, the closed form of the z-transform is 

u / \ _ z (fl! - £ 2 e ) + 

2 ' f ! (e^ ,T )- 2 -'p T ( e^) t . 


(63) 


namely. 


It is observed that H-j(z) and H 2 (z) are of the same form. 


H(z) = 


A 0 f Aiz 


-1 


1 + B-j z” 1 + B 2 z" 2 


(64) 


This transform can be realized in several ways Assume the input 
sequence is x(n) and the output sequence is y(n), as shown in Figure 
36. In the z-transform domain 


y(z) 


H(z)X(z) = 


A 0 + A 1 z 


-1 


1 + B-j z' 


„ -- 2 xC2) 


+ B 2 z 


(65) 
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x(n) 

H(z) 

y(n) . 


, - -- ” r 


Figure 36. General Block Diagram of Digital Simulation 

for Acceleration 

where, for causal filters, z -1 represents a time delay of T seconds, 

T being the computer update time. In the time domain, 

y(n) = - (b-j y(n-l ) + B 2 y(n-2)) + A Q x(n) + A ] x(n-l) . (66) 

Three different realizations of this type of transfer function 
are shown in Figures 37 through 39. In Figure 37, the direct recursive 
realization is shown. The constants are those of H(z) in (64). The 
Canonic Realization is shown in Figure 38, where the constants are also 
those in (64). In both of these realizations, there is a severe accu- 
racy limitation on the constants {A^ and {B^. As an alternative, 
the Parallel Canonic Realization is shown in Figure 39. This reali- 
zation implements the representation 

§1 

H(z) = L - TT + =r , (67) 

i + &, z i + e 2 z 

which is a partial fraction expansion of H(z) in (64). This implemen- 
tation is recommended since the accuracy on {<$•} and in the Parallel 
Canonic Realization is not nearly as severe as the accuracies in the 
other two candidate realizations. 

In this section, the similarities of the digital algorithms for 
simulating the effects of acceleration orthogonal to and along the LOS 
have been presented. The closed loop transfer functions between accel- 
eration and the corresponding tracking error are given for all tracking 
loops in the Ku-band radar. Finally, specific implementation realiza- 
tions are considered, and the recommended one is specified. 




Figure 38. Canonic Realization 
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Figure 39. Parallel Canonic Realization 
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Range Rate Logan thmic Discriminant Error Simulation 


In order to simulate the range rate logarithmic discriminant 
error, it is important to examine the nature of the problem. The range 
rate is estimated during track by using a set of 32 doppler filters 
and a log discriminant generator. The mean of the output of the log 
discriminant generator [E(z)] was found in [1, Appendix G] to be: 



Figure 40 illustrates the variation of E(z) as a function of s with SNR 
as a parameter. Since the SNR is not known, a given value of E(z) does 
not result in a single value for s because the mapping is not one-to-one, 
which gives rise to the range rate logarithmic discriminant error. In 
any implementation, a given one-to-one mapping has to be established. 

The error can then be found as a function of the signal-to-noise ratio. 
The maximum error found in the expected range of operating signal-to- 
noise ratios can then be used as a bound for this error. Two methods 
will be presented here. The first is being implemented by Hughes Air- 
craft Corporation, and the second is included as an alternative. The 
simulation has two parts: The first is to find the value of af d given 
a certain range rate v^, and the second is to find the error using a 
given method of implementation. 
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Normalized Doppler Offset, c 

Figure 40. The Mean of the Logarithmic Discriminant Output Versus 

Normalized Doppler Offset 
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4.3 1 Finding the Doppler Offset Afj 

Since the doppler filters remain fixed in the tracking mode with 
a spacing of (Bp/2) Hz apart, the doppler reading is obtained using the 
filter whose center frequency is closest to the actual doppler. The 
two adjacent filter readings are used as inputs to the logarithmic dis- 
criminant. Figure 41 illustrates the location of the 32 filters. The 
number between two adjacent dotted lines represents the filter which is 
used for doppler measurement when the actual doppler lies between the 
two lines. To obtain the doppler offset Af^, we assume that the velocity 
(range rate) R is positive; the actual doppler is then equal to 


f 


d 



(69) 


where f is the RF frequency and c is the velocity of light in free 

^ t 

space. For R in ft/sec. 


f = 2 (13.885) 10 ^ = 28.2 R (Hz) 
a 9.83569 x 10° 


We select the filter n to be such that 


okxotnalpaseb 

OF POOR QUALU-* 


(2n-l ) Bp 
2 




( 2n+l ) Bp 

J 

2 


n = 1 ,32 . 


(70) 


If n.[ is the particular n satisfying (70), then 


Af 


n l b f 


- f. 




(71) 


This value of Af^ is used to calculate ? = Af^/Bp There are several 
points to observe: 

(1) ? is a positive number which is less than 1/4. 

(2) The absolute value used in (71) is due to the symmetry 
of the problem. 

(3) For negative range rates, the absolute value of the 
velocity is used in the calculations resulting in the number n-j. 

Since the filters are arranged in a circular configuration, the filter 
(32-n-|) is used instead of n^ to obtain Af d in (71). 

(4) The filter numbers 32 and 0 represent the same filter. 
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4.3.2 Methods for Finding the Error 

Assuming that a one-to-one mapping between E{z) and ? exists, 
then for a given reading E{z), a value e-| can be found using that 
mapping. The range rate logarithmic discriminant error then becomes 

A \d = (?- s-j) Bp/28.2 ft/sec. (72) 

Finding ^ depends on the particular mapping. Two such mappings are 
considered. 

(1 ) HAC Mapping 

This is the actual mapping being implemented by Hughes. 

A straight line with a slope equal to that of E(z) at £ = 0 with SNR=® 
is used to obtain c-j- This is the straight line mapping as shown in 
Figure 40 The slope (m) of the line can be found from (68) to be 

m * -J- cot (i^r) • ( 73 ) 


The resulting range rate logarithmic discriminant error is shown in 
Figure 42 along with HAC values for comparison purposes. 

(2) Infinite Signal-to-Noise Ratio Mapping 

Since the ranges of primary interest are close ranges, 
E(z) at SNR = “ can be used as a reference (mapping). The value of e-j 
can be found from (68) by equating SNR to infinity, resulting in 


4 


1 


1 

2 


2e E(z)/2 


- 1 . 


(M) 


The actual error at a given SNR can be found using (72). More 
accurate methods can be suggested to minimize the range rate log dis- 
criminant error. All of these methods, however, require major design 
modifications, which are neither warranted nor recommended at this 
stage. 





Range Rate Log Discriminant Error (ft/sec) 
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Figure 42. Range Rate Log Discriminant Error, Using HAC Method 
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4.4 Simulation of Target Effects 

Models are given in this section to take into account the fact 
that a true radar target is not a point target. Various models have 
been prepared in the past, but the simplest and possibly the most effec- 
tive are those described by Barton [11-12]. Variations of those models 
are described and recommended for the Ku-band radar simulation. The 
target effects on each of the tracking loops are considered individually. 
In each case, an analog computer model is given which involves broad- 
band noise being passed through a first-order filter. The model using 
a digital recursive algorithm is presented in Appendix A of this report 
and applies directly to each case presented here. 

4.4.1 Target Effects on Range Errors 

A block diagram of the analog model which represents target 
effects on range errors is given in Figure 43. The output ARy(t) is 
the range error induced by the target, causing time scintillation of 
the transmitted pulse. The total power or, equivalently, the standard 
deviation of this error is approximated as 

° 4 R t ■ a r L r (75 > 

where the units of are the same as L r (feet or meters) and a r is 
a proportionality constant in the range 

0.2 < a r < 0.5 . (76) 

This says that the standard deviation of range errors due to target 
effects is 20% to 50% of L^, where L r is the effective target dimen- 
sion along the LOS. This is a value which must be approximated and 
can also be time varying. 


n r (t) 



1 


T r s +1 


AR t ( t ) 

— ► 


Figure 43 Analog Model of Target Effects on Range Errors 
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A simple and satisfactory model of the target (in a plane) is 
a rectangle with minimum dimension L n and maximum dimension L . 

min niaA 

We assume the target is rotating in the plane created by the LOS by 
an amount radians/sec. Then the effective target distances along 


and perpendicular to the LOS are 


L r (t) “ Lq [1 + y cos (» p t)] 

(77a) 

and 


L„(t) = L n [1 + y sin (to t)] , 

(77b) 


respectively, where 


and 

Y 


. L + L * 

A max mi n 


* L „ax - Spin 

L + L . 
max min 


(78a) 


(78b) 


In this model (shown in Figure 44), it is assumed that the rota- 
tion is small compared to the reciprocal of the tracking loop bandwidths. 
This is certainly the case with any target for which there is expected 
rendezvous and docking. As a result, L (t) and L (t) need not be updated 
every T sec (computer update time) and could be updated every 10 T to 
20 T. 


The time constant of the target along the LOS in Figure 43 
is given by 


T 


r 


4u L 5 
r r 


(79) 


where X is the wavelength. If the target's relative rotation with 
respect to the LOS is identically zero, then this time constant goes 
to infinity, and the range error due to target effects is a constant 
bias error. For this time constant and the standard deviation in (75), 
the one-sided power spectral density (PSD), W^., of the broadband noise 
n f (t) is given by 
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3 L A 2 

= — - (units of distance) sec . 


( 80 ) 


A block diagram of the computations necessary for this simula- 
tion of target effects in range is shown in Figure 45 The stochastic 
differential equation for ARj.(t) is 

AR T (t) = (^-) AR T (t> + (J-) n r (t) . (81) 


To convert this to a discrete time digital simulation, follow the pro- 
cedure at the end -of Appendix A of this report. In Figure 45, the 
factor l/v'IT accounts for the improvement on the target effect errors 
which results from frequency diversity. 


4.4.2 Target Effects on Range Rate Errors 

A block diagram of the model of target effects on range rate 
errors is shown in Figure 46. The standard deviation of the error is 
approximated by 


2a v L x“r 


Hz , 


(82) 


where a y is a proportionality constant which, for most targets, is 
within the interval 


0.2 < a y < 0.4 . (83) 

In terms of range rate, the standard deviation is 

°v “ a v*"x u r ( units of distance)/sec . (84) 


The time constant of the target effects on range rate, t y , is 
approximately equal to 


T 


V 


4 a) L 
r x 


sec , 


(85) 


where uy is the same rotation rate of the previous section m rad/sec. 
If Up = 0, there are no target effects on range rate. 




Figure 45. Simulation of Target Effects in Range 
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n (t) 

v - fc. 

1 

AR (t) 

► 

v + 1 

W 


Figure 46. Analog Model of Target Effects on Range Rate Errors 


For the standard deviation given in (84) and the time constant 
in (85), the one-sided PSD of the noise n y (t) is 

2 2 

W v = a y L x *“ r (units of distance) /sec. (86) 

The block diagram which depicts the computation necessary to 
simulate the target effects on range rate is given in Figure 47. The 
effect of frequency diversity (a factor of 5 improvement in variance) 
is also shown. 

The stochastic differential equation for AR^(t) is 

■ (17) + §) "v (t > ■ < 87 > 

which can be simulated digitally via Appendix A. 

An example calculation of the standard deviation is: 

Example : a y = 0.35 

L = 30 ft 

A 

tii = 2 deg/sec 

„ . ( = 0.37 ft/sec (without frequency diversity) 

U £R ) 

T( = 0.17 ft/sec (with frequency diversity) 

Note that the target effects on both range and range rate are 
independent of the true range. It will be seen that this is not the 
case with the angle and angle rate estimates. 

4.4.3 Target Effects on the Angle Tracking Loop 

The target effects on the angle tracking loop are similar to 
those of range and range rate presented above with the exception that 
the effect depends on the range to the target. The block diagram is 
shown in Figure 48. The standard deviation of the angle error due to 
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target effects is approximated by 

0 A0 = “T“ rad1anS ’ 


( 88 ) 


where the proportionally constant a 0 is also in the range 

0.2 < a < 0.4 (89) 

u 


and L„ and R are in the same units. The time constant x a of the analog 
x y 

model is equal to 


T 


e 


7 p- sec , 

4 “a L x 


(90) 


where to (in rad/sec) is angular rotation in the plane perpendicular 

a 

to the LOS. The corresponding PSD, W Q , of the broadband noise n Q (t) is 


W 


0 


2 

a 9 L X X 
R 2 to 

a 


sec . 


(91) 


The necessary computations are shown in Figure 48. The stochastic 
differential equation for A0(t) is 

. M) „(t) ♦ (i) n 9 (t) , (92) 

9 9 

which is simulated, as are the others, via Appendix A. 

The standard deviation in (88) is plotted in Figure 49 for cases 
that are expected on the Orbiter. It is seen that the target effects 
can be expected to cause the angle error to be over the specification 
value at relatively large ranges For the values chosen, the range at 

which the antenna field of view is 10% and 100% filled is indicated on 
Figure 49 using the computation described in Figure 50. 

The results indicate that reliable angle tracking information 
cannot be expected at ranges above several thousand feet when L = 60 ft. 

X 

A rule of thumb is that, when the antenna beam is filled by more than 
10%, angle tracking difficulties should be expected. In the case of the 
Orbiter, however, this percentage can be increased somewhat because of 
the slow rotation rates that are anticipated. The solution to this prob- 
lem is described in the introduction of this report. 








Field of View, d (m) 



Figure 50. Range at Which the Antenna Field of View is Filled 
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4.4.4 Target Effects on Angle Rate Estimation 

The angle rate error caused by target effects is best estimated 
from the development in [1, Appendix D], wherein the rates of the 
standard deviations are given by 

cr* 2irf , 

-T - n (sec)" 1 . (93) 

where f = natural frequency of the angle tracking loop 
C = damping factor of the angle tracking loop . 

This can also be written as 

. 2 * f n 

Ji 2s) Z (»f n ) 2 +1 

where t is the time constant of the angle tracking loop. For R<1.9 nmi, 
the present system parameters [2] are 

2t = 5.4 sec 

f n = 0 12 Hz, (95) 

for which the term in brackets in (94) is equal to 0.332 (sec)”^. 

Since the angle rate estimate is derived directly from the angle 
tracking loop (see [1, Appendix D] or section 4.1 of this report), the 
target effects for angle rate will have the same power and time con- 
stants as employed for the angle tracking loop. Angle rate errors can 
be determined directly from angle errors as shown in Figure 51. 

Sample computations of the angle rate error are shown in Figure 
52, where it is noted that the specification value on standard deviation 
is not attained below several thousand feet for most target sizes. 

An alternative solution to this difficulty has been proposed and 
is described in the introduction of this report. 

4.5 Overall Description of Ku-Band Radar Simulation 

In the previous parts of this section, the necessary components 
of the Ku-band radar simulation were presented. In this section, the 
components are assembled and certain additional errors are taken into 
account. 


(94) 
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Figure 51. Angle Rate Errors From Angle Errors Due to Target Effects 
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Figure 52. Angle Rate Errors Due to Target Effects 
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4.5.1 Range Error Simulation 

The range error simulation block diagram is given in Figure 53. 

The target effects, thermal noise and amplitude scintillation errors, 
and acceleration errors have been previously diagrammed. The component 
bias is a range error which results from time delay variations of the 
RF signal through the various components in the front end of the radar. 
These variations are due to temperature, aging, and the like. The 
maximum value of these components [2] has been estimated as 

L = | ±28' | max for R< 2500' 

L = |±77* i max for R> 2500' . (96) 

The variation is very slow. Therefore, we propose to model the component 
bias as a random variable, either Gaussian or uniform. In the Gaussian 
case, set 3o=L, with L given in (96). In the uniform case, the proba- 
bility density function extends from (-L,L), and the variance is 
o= L//3". 

A second additional error component of range error is due to 
data staleness. Using the update times of the computer involved with 
the Ku-band radar, a stochastic model can be developed similar to those 
used in the target effects simulation. A simpler alternative is to 
periodically use a sample of a Gaussian noise generator whose 3o value 
is equal to 1/2 the maximum data staleness value and whose mean has 
the same value. This is depicted in Figure 54, where M corresponds to 
the maximum data staleness time. 

A sample computation of range errors in the proximity operations 
region is shown in Table 9. It is seen that, at short ranges, the com- 
ponent bias and target effects constitute the major part of the range 
error. The acceleration effects are small, and the effects of data 
staleness were not taken into account in the computations of Table 9. 

45.2 Range Rate Error Simulation 

A general block diagram of the range rate error simulation is 
shown in Figure 55. All contributions to the range rate error have 
been presented in the previous sections A sample computation is shown 
in Table 10. The maximum range rate error due to the logarithmic dis- 
criminant error would occur only at lower signal-to-noise ratios. At 
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Figure 53. Range Error Simulation 
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Figure 54. PDF for Data Staleness Model 


Table 9. Typical Case of Proximity Operations Range Errors 


1 . 

Thermal Noise 

a * 0.1 ft 

2 . 

Target Effects 

a = a L ft 
r r 

3. 

Component Bias 

cr = 28/3 ft (Gaussian) 
= 28//3 ft (Uniform) 

4. 

Error from Acceleration R 

Small 

5. 

Data Staleness 

Gaussian Random Variable 

Example: a r = 0.3 ® total 

L r = 30 ft 

= 18.4 ft (Gaussian 
= 25.3 ft (Uniform) 

Table 10. Proximity Operations Range Rate Errors 

1 . 

Thermal Noise 

a < 0.1 ft/sec 

2. 

Target Effects (Example in 
Section 4 5.2} 

a ~ 0.17 ft/sec 

3. 

Ln Discriminant Error 

Max 0.4 ft/sec 

4. 

Acceleration Effects 

Negligible 


“total ” 0 ' 5 ft/sec 




Figure 55. Range Rate Error Simulation 








96 


close ranges, the SNR will be quite high, so that the logarithmic dis- 
criminant error will be much less than the maximum indicated in Table 
10. In the absence of significant target rotation, the Ku-band radar 
can be expected to provide satisfactory range rate data down to a few 
hundred feet. 

4.5.3 Angie and Angle Rate Error Simulation 

The general block diagram of the angle and angle rate error simu- 
lation appears in Figure 56. All contributions to these errors have 
been discussed in previous sections. The primary contribution at close 
ranges is the target effects, plots of which are presented in section 
4.4.4. 

4.5.4 Additional Comments 

The acceleration effects and thermal noise effects took into 
account the transfer function of the Ku-band radar. When the target 
effects were taken into account, the effects of the transfer function 
were not included. This would need to be done if the power spectral 
density of the target effects were wider than those of the radar track- 
ing loop. Since the target rotation is expected to be very small, the 
bandwidth of the signal will not be more than 0.1 Hz or -less. There- 
fore, target effect errors will be direct errors, and radar filtering 
will not alter them. 

As a final remark, it should be noted that no range gates are 
used at close range ( R < 0 . 42 nmi). Instead of the straddling loss, 
there is a substantial signal loss due to the two-pole Butterworth 
filter when the range is short and the pulse width is 0.122 ysec. The 
response to such a pulse is shown in Figure 57. Since only one sample 
per pulse is taken at short ranges, the maximum possible power loss is 
found to be 6.5 dB. This loss should be included in the radar equation 
when computing the signal-to-noise ratio which is to be used in the 
simulation computations. 
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Figure 56. Angle and Angle 
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5.0 ELECTRO-OPTICAL DOCKING TECHNIQUES 

This section describes the results of a preliminary tradeoff 
analysis conducted to identify viable electro-optical (EO) docking 
techniques to supplement the Ku-band radar at short ranges from 
approximately 1000 ft down to 25 ft. 

Section 5.1 presents an analysis of the geometrical relations 
involved in the co-planar case, including an error-sensitivity analysis, 
while Section 5.2 contains a description of the various EO system con- 
cepts considered for both cooperative and noncooperative satellites. 

The primary device for implementation of these techniques is the 
laser. Further study and analysis of the proposed methods are needed 
to confirm the tentative conclusions given here. 

5.1 Analysis of Ranging Techniques 

In this section, the geometrical relationships used to compute 
the "perpendicular range" from the baseline of the docking station 
(Orbiter/Shuttle) to the satellite (target) are developed. Two basic 
methods are defined- (1) tri angulation (angle measurement) and (2) 
range-only measurements. Possible combinations of these are also 
noted. Various electro-optical concepts utilizing these techniques 
are described in Section 5.2. 

5.1.1 Tnangulation 

The tnangulation method requires two angle trackers situated 
at the extreme ends of the Shuttle baseline (to maximize range measure- 
ment accuracy). If the distance between them is accurately known and 
the appropriate line-of-sight (LOS) angles to the satellite are measured 
by the trackers, the perpendicular distance from the baseline to the 
satellite can be determined by trigonometry When the tracker angles 
are measured in a coordinate system that is fixed with respect to the 
Shuttle, the "range triangle" will not, in general, lie in the reference 
plane Although general, three-dimensional relationships would be 
employed in the real docking system (i.e., they would be programmed 
into the docking computer), only the much simpler co-planar case will 
be considered here, since it suffices for preliminary analysis. 
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The co-planar case is illustrated in Figure 58, where a and b 
are legs of the triangle, a is the base (baseline), and the angles 
and $2 are measured from the normal to the baseline (i e., from the 
y-axis) as shown. The dashed line R is the perpendicular distance to 
the satellite and uniquely defines its range from the station. 

By straightforward trigonometry, it can be shown that 

n cos 4>-| cos $2 
R = sirTffj +f 2 ) • 


which involves the two variables <f>-| and <j> 2 in addition to the constant 
£. To avoid unnecessary complexity, we restrict further consideration 
to the special case when the satellite is equidistant from the two 
trackers and approaches along the perpendicular R. Then, <f>-j = = <j>> 

and 

2 2 
p - k cos & - k cos ^ 

sin 2<p 2 sin (f> cos <j> 


R($) = D cot 4> , 


(97) 


where D=a/ 2. Equation (97) is seen to be correct by inspection of 
Figure 59. The notation R(<{>) signifies that R is a function of the 
measured variable 4>. 


Error Analysis for R 

We are interested in the accuracy of the R determination in terms 
of the uncertainties in the measurements of <j> and D. To this end, con- 
sider the variations 6R^ and 6R^ resulting from independent variations 
in (j> and D (or a), respectively. Thus, from (97) and Figure 59, 




D esc 4» 


2 

R +D 
D 



D <5 <j> 

2 

sin <j> 


6Rg = cot <f> <5D 
6R D = 2D 6£ ‘ 


and 
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Figure 58. Co-Planar Case of Tri angulation 



Figure 59. Special Case of Co-Planar Triangulation (b = a) 
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If and 6D are regarded as independent random variables, then 
the rms variation in R is 

5R ms ‘ ./< 5R / + (sR / 

6R rms = ^ y (R 2 +D 2 ) 2 (6<J>) 2 + R 2 (6£/2) 2 . (98) 

Hence, it follows that 

J d 2 (sr) 2 - r2 ( 6D ) 2 , , 

S4> = r A S . (99) 

Thus, 6<{>(rad) is the angle accuracy required to achieve a range pre- 
cision (<5R) = 1 ft at the range R. 

rms 

Equation (99) is plotted as a function of R in Figure 60 for 
D=30 ft (a= 60 ft), = 0.1 ft, and 6R rms =0.5, 1, and 2 ft. These 
values are typical of those anticipated on the Orbiter. Since the 
range precision (<5R rms ) is held fixed, the allowable angle error 6<j> 
decreases as the range increases and goes to zero at a range R m , 
where from (99), 

A ^ SR Vms 
6 D 


When R>R fn , the range error exceeds ( fiR ) rms even for 5<j>=0. 

In practice, however, the allowable range error increases with 
the range and has its minimum value at minimum range. Taking the mini- 
mum range to be 25 ft. Figure 60 shows that the required angular preci- 
sion to ensure (<SR) rms values of 0.5, 1, and 2 ft at this distance is 
9.8, 19 7, and 39 3 milliradians, respectively All of these values 
are achievable with closed-loop trackers, and at least the larger two 
are also compatible with manual trackers. If these values are held 

fixed (as in a real system), then the range error increases with R as 
shown in Figure 61. Although the error is quite large beyond 100 ft 
or so, it should be tolerable at such distances, provided the satellite 
is not closing too rapidly. Also, these values can be reduced by 




Figure 60. An 
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further restricting the angle error tolerance 6<t>, particularly if closed- 
loop trackers are employed. 

Error Analysis for R 

The corresponding formulas for range rate R and range rate error 
6R follow in similar fashion after taking the time derivative of (97): 

R = D cot $ 

R(<M) = “0 esc <f> <j> ( D — 0) . (100) 


In this case, variations in R can arise from variations in the measure- 
ment of <}>, <f>, or D. Hence, from (100) and Figure 59: 

( 6R) = - 2D esc 4> (- cot <j> esc <{>) $ 6<f> = - 2(- D esc 4 <j> $) cot <j> 6$ 

= - 2 R( R/D) 6<j> 


(«R)j 


= - D CSC <t> 6(f) 


so 


and 


r2 +D 2 
— D — ^ 


(6R) d = - esc $ <5D 

= R(«D/D) , 

(«R) ras - J(s R) t 2 + t«R);+ UR) 


|j2R/D) 2 (6<f.) 2 + ( 6D/D )®J A 2 + ( -- p — ) (4) 2 (101) 


6<j) “ 


D 


^ — JJ^ms ' [(2R/ D ) Z (^) Z + (WD) 2 ]R 2 
R + D K 


( 102 ) 


5.2 Direct Rangefinder (No Angle Measurement) 

This method employs two independent rangefinders to measure the 
legs a and b of the range triangle. If the baseline & is also known. 
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the perpendicular range to the satellite can be determined. Although 
no angle measurements are used, an automatic system would require closed- 
loop angle trackers to point the lasers continuously at the satellite. 
These trackers could also provide an independent and redundant range 
determination based on the triangular method (section 5.1). However, 
if the lasers are pointed manually, no trackers are needed. 

The desired expression for range R in terms of a, b, and z is 
easily deduced (in the co-planar case) from the geometry of Figure 58. 

By the law of cosines, 

cos e 


a 2,_2 k 2 

- a -f £ - b 
2 a i 


and 


2 / „i 2 

- a - (a cos 6 ) 



003) 


In the special case when 0 = 6 , a=b, 

r - ^ - {0 - • < ,o4a > 

where D = z/Z. 

Error Analysis 
Squaring (104a) gives 

R 2 = a 2 -D 2 . (104b) 

Hence, 

R ( <5R) .j = a 6 a 
R( 6 R) D = -D 6 D, 


( 6 R) 


rms 


so 


Jis R a ) 2 + (SR D ) 2 = lJa Z {6a) 2 + D 2 ( 6 D ) 2 


(105) 
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and 

«a = D ^ 6D > 2 • 


(106) 


where 6a is the rangefinder precision required to ensure a perpendicular 
range accuracy of (SR) ft at range R, for specified 6D. 

I 

The angle error can be similarly found to be equal to 


(63) 


rms 


J(6D) 2 + (D/a) 2 (6a) 2 
R 


If D is known without error, then the rms angle error can be written 


as 


(63) 


D 6a 


rms 




rjr 2 +d 2 


5.3 Description of Electro-Optical Concepts 

In this section, various electro-optical concepts for implement- 
ing the ranging techniques of section 5.2 are described, and recommended 
concepts are selected for cooperative and noncooperative satellites. 

All of the concepts considered derive the necessary geometrical infor- 
mation by measuring the distance or angular position of a fixed refer- 
ence point on the approaching satellite by means of two E0 sensors 
positioned at either the same or opposite ends of the platform baseline. 
The satellite reference point may be either self-luminous (beacon con- 
cept) or illuminated by a laser (or lasers) located on the platform 
(laser/retro or laser spot tracker concepts). 

Purely passive concepts such as tracking the satellite by solar 
illumination or thermal (black-body) radiation have not been examined 
in detail because the former is too variable and unreliable while the 
latter presents a radiation signature that is too weak and spatially 
diffuse for accurate terminal guidance, particularly in the presence 
of strong celestial or earth backgrounds 

The sensors considered are all designed to operate over a 
limited field of view. While some are nulling devices, some means 
must be employed to steer the sensors continuously This is because 
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the angular position of the reference point changes during satellite 
approach. The same capability is required for the lasers, which pro- 
duce an inherently narrow beam that must be kept continuously on the 
illuminated reference point in order to ensure an accurate range or 
angle measurement. 

The necessary steering capability can be provided by one or both 
of the following methods: (1) manual steering of either one or two 

television systems fitted with viewing reticles and electrically slaved 
to the sensors and/or lasers, and (2) automatic steering by means of 
closed-loop trackers. The first method is essentially an open-loop 
approach which is limited by the achievable tracking precision of a 
human operator. Since suitable TVs are already available, this is an 
attractive possibility, provided the tracking jitter is not excessive. 

The second method can provide much better tracking precision 
than manual tracking (a particularly important advantage with the tri- 
angulation technique). However, it is more complex since it requires 
a complete closed-loop servo control system, and it does not provide 
automatic acquisition unless a special acquisition mode is designed 
into it. (A hybrid approach, combining manual TV acquisition with 
handover to closed-loop tracking, is also a very practical solution 
to this particular problem.) 

Several types of two-axis tracking sensors can be considered for 
closed-loop trackers. Among these are silicon quadrant or continuous- 
position sensing detectors (for visible and near-IR wavelengths) and 
moving-reticle type star trackers using a photomultiplier tube. All of 
these sense the displacement of an optical image of a point source rela- 
tive to two orthogonal axes (e.g., azimuth and elevation) and thereby 
measure the angular position of the satellite reference point with 
respect to the instantaneous line of sight (LOS). The resulting terror 

i 

signal can be used to command a closed-loop servo system. Shaft angle 
encoders would be needed to measure the net LOS rotation relative to 
coordinate axes fixed in the docking platform, as required for the tri- 
angulation technique 

Although closed-loop TV tracking is technically feasible, it has 
been excluded from detailed consideration because of its complexity and, 
cost relative to other methods. It should be noted, however, that TV 
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trackers are especially vulnerable to image blooming and tube failure 
resulting from excessive Image brightness such as might occur due to a 
strong sunglint off the satellite. 

5.3.1 Cooperative Satellites 

A cooperative satellite is defined to be one that carries with 
it one of the active or passive components of a complete docking system. 
For the E0 system considered here, this component may be either a beacon 
or a retro-reflector array. The beacon may be activated either from 
on-board the satellite or remotely, by a signal transmitted from the 
docking platform. {In the latter case, it would be regarded as a 
transponder.) Pulsed operation is preferable to a CW mode, both to 
conserve average power and to aid in background discrimination. If a 
retro array is used, a suitable laser or other high intensity illumina- 
tion source must be made available on the platform. 

Both the beacon and the reflector array provide a tracking/ 
ranging reference point that has a fixed location on the satellite, 
contributing to more stable tracking than can be achieved by tracking 
a laser- illuminated spot on the skin of the satellite. This tracking 
reference is also generally much brighter than the illuminated skin and 
therefore produces higher signal -to-noise ratios in the tracker. 

Beacon 

The beacon concept utilizes ranging by the triangulation method. 

It was shown in section 5.1 that this method is capable of achieving 
sufficient range precision at minimum range, using closed-loop trackers 
of ordinary accuracy. (The adequacy of manual TV tracking is less 
certain, but this technique may also be viable.) If the precision 
achievable beyond the minimum range is also deemed acceptable (see 
Figure 61), then the beacon is the preferred concept for cooperative 
satellites because (1) no active sources are needed on the docking 
station; (2) manual acquisition via TV is easily implemented by slaving 
the tracking sensors to the TV line-of-sight (the sensor's field of 
view must be large enough to accommodate LOS differences due to parallax) 
(3) problems associated with pointing of a laser illuminator or range- 
finder are avoided; and (4) the beacon provides a strong, reliable signal 
exhibiting less fluctuation than is present with other methods. 
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The primary drawback of a beacon for a cooperative satellite 
is that it consumes some spacecraft power. It also requires timely 
maintenance to ensure that it operates reliably at the time of docking. 

The following two types of beacons are viable candidates: 

(a) Xenon Aircraft Beacon (Recommended Concept) . The Xe beacon 
is recommended over other alternatives subject to the proviso that 
excessive spectral filtering is not needed to contend with sunglint 
background. If such filtering should be found necessary (through later 
analysis) and results in unacceptable signal loss, then the Gallium 
Arsenide beacon described below would be favored instead, because of 
its intrinsically narrowband output. 

Aside from the question of the filtering required, if any, the 
Xenon beacon is a highly attractive concept for cooperative satellites 
because of its wide beam, visibility, TV-compatible output spectrum, 
and because it is a standard catalog item requiring only minimal inter- 
facing design. A small, pulsed beacon drawing a few watts of average 
power should be sufficient for docking control out to a range of at 
least 1000 ft. It may be advisable to include a simple reflector to 
confine the beam to the forward hemisphere, thereby increasing the 
intensity relative to the isotropic mode. 

(b) Gallium Arsenide (GaAs) Beacon . This device consists of a 
close-packed array of GaAs laser diodes with an associated pulse and 
simple optics to project the beam. Arrays and pulsers are commercially 
available Uncoded GaAs diodes emit at a wavelength of about 0.9 m 
(900 nm) in the near-infrared, unless special doping is used to shift 
the output to a shorter wavelength. Aluminum-doped GaAs (GaAlAs) diodes 
emitting at wavelengths of 0.8 pm or less can be obtained from at least 
one commercial vendor at the cost of somewhat reduced peak power per 
diode. This expedient may be necessary to make the beacon visible to 
the TV for tracking and/or acquisition. Even so, a TV with an extended- 
red response may be required. 

The raw output beam from a GaAs diode array tends toward a conical 
shape having a full -angle divergence on the order of 30 degrees at the 
50% intensity points. By the use of a simple optical mixer or "scrambler," 
the output can be made highly uniform, and it is then particularly well 
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suited for a wide-angle beacon, with little or no additional optics. 

The requisite intensity is achieved by employing a sufficient number 
of diodes in the array. Since the individual diodes are typically only 
a few mils wide in the longer face dimension, the overall size of the 
array is quite small, and it is relatively lightweight. 

Cooled GaAs arrays offer substantially higher average power and 
efficiency and also emit at a slightly shorter wavelength. However, 
for the present application, these advantages are more than offset by 
the additional complexity and maintenance requirements of the associated 
cryogenics. Hence, an uncooled array is recommended unless further 
analysis shows that adequate intensity cannot be achieved without 
cooling.* 

In comparison to the Xenon beacon, the GaAs beacon offers the 
advantage of an intrinsically narrow output spectrum (helpful for back- 
ground discrimination), compact design (if cryogenic cooling is not 
required), potentially lower power consumption (depending on the array 
size and beamwidth employed), and possibly longer service-free intervals. 

On the other hand, the GaAs beacon is power-limited to a relatively 
narrow beam and hence imposes some constraints on satellite orientation, 
and its infrared output is poorly matched to the spectr al re sponse of 
the TV. Further analysis will be needed to make a definitive choice 
between these alternatives. 

Laser/Retro 

This concept utilizes a retro-reflector array on the satellite 
and one or two lasers on the docking platform. Several types of reflec- 
tors may be considered. If the array is comprised of individual corner- 
cube reflectors, the return beam will be highly directional and will 
only be visible to a sensor co-located with the laser. Other devices 
for enhancing the laser return include glass-bead type reflectors and 
reflective tape, which tend to be somewhat less directional than corner- 
cube reflectors. A fairly broad return beam is necessary if a tracking 
sensor is to be located at the opposite end of the baseline from the laser. 


* # 

Radiative cooling to space may be an effective alternative to 
cryogenic cooling and should be considered. 
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Because of their narrow beams, the lasers must be pointed con- 
tinuously at the retro as the satellite closes in range. Steering can 
be accomplished by either of two methods. First, the beam divergence 
may be increased sufficiently {using external optics) to permit manual 
TV steering, provided that a sufficient SNR can be maintained with the 
broadened beam. Alternately, a closed-loop tracker (or trackers) can 
be employed. In this case, the beam(s) must be nutated or conically 
scanned about its axis to generate a tracking signal, while its diver- 
gence must be compatible with the tracking accuracy, jitter and lag of 
the tracker, to ensure that the beam does not drift off the retro. 

With either open-loop or closed-loop tracking, acquisition can 
be achieved by manually pointing with a wide beam, using radar data to 
define the uncertainty volume. Acquisition will be signaled by the 
appearance of a bright spot on the TV monitor. The beam divergence can 
then be reduced, if desired, to improve the SNR for closed-loop tracking. 

The laser/retro concept can be implemented in several ways, the 
principal ones being the following: 

(a) One laser illuminator/ rangefinder, one angle tracker . In 
this method, a CW-modulated laser serves as both a rangefinder trans- 
mitter, with a co-located receiver, and as an illuminator for an angle 
tracker, also co-located with the transmitter. This hybrid system 
measures the length and angle of one side of the range triangle and 
hence represents a combination of the rangefinding and tri angulation 
techniques. 

(b) Two laser rangefinders . This range-only approach measures 
the lengths of the two legs of the triangle, as described in section 5.2. 
If corner-cube retros are used, the two return beams will be isolated 
from each other in angle (each laser/sensor "sees" a different retro- 
reflector of the array) and lasers of the same wavelength can be employed. 
Otherwise, different wavelengths and spectral filters should be used 

to prevent possible mutual interference from overlapping beams. By 
the addition of one or two angle trackers, as in (a), the system can 
provide a redundant or backup determination of range by triangulation. 

(c) Two laser illuminators/trackers . This variation uses the 
lasers only as illuminators for their associated angle trackers, and 
perpendicular range is obtained by triangulation. Two illuminators 
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are needed only if corner-cube or other highly directional retro reflec- 
tors are used on the satellite. 

The laser /retro concept shares most of the virtues of the beacon, 
with the additional advantage that the satellite is required to carry 
only a passive device that consumes no power. Of course, this simply 
transfers the active device to the docking platform. Another signifi- 
cant advantage is that an illuminator or rangefinder laser can also be 
used for tracking/ranging noncooperative satellites, as described in 
the following section, whereas beacon systems cannot. 

The major disadvantages of the laser/retro concept are the poten- 
tial difficulty in acquiring the retro with a laser and the possible 
requirement for a closed-loop laser pointing and tracking system, which 
would substantially increase the complexity and cost relative to the 
beacon system. 

Search! ight/Retro 

This concept envisions high-intensity, noncoherent, broad-beam 
sources on the platform to floodlight the entire satellite as it 
approaches. Little or no steering of the sources would be required. 

The retro-reflectors would appear as relatively bright points against 
a lesser background of random glints and diffuse reflection from the 
skin of the satellite. By properly thresholding the tracking sensors, 
it should be possible to discriminate against this background In other 
respects, the concept is identical to the two-laser llluminator/tracker 
method described above. 

Possible sources include the Xenon and GaAs beacons discussed 
earlier, situated on the platform rather than on the satellite. The 
major uncertainty concerns the source power necessary to provide a 
sufficient retro-return at maximum range. Otherwise, the concept is 
an attractive one, since steering requirements are minimal. 

5.3.2 Noncooperative Satellites 

As noted earlier, passive tracking methods based on solar reflec- 
tion or thermal emission from the satellite do not appear to be viable 
candidates for this application. Consequently, only active EO concepts 
using lasers (or possibly other high intensity sources) on the docking 
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platform need to be considered for noncooperative satellites. Further- 
more, the requirement for a well-defined tracking point excludes broad- 
beam sources that would floodlight the satellite, since the multiple 
glints produced would only confuse a tracking sensor. Thus, the only 
viable approaches for passive satellites are those involving lasers 
on the station, since they are capable of illuminating a small spot at 
ranges of 1000 ft or more. 

The options available with lasers are basically similar to those 
for the laser/retro concept, as described above, except that the beams 
are reflected from the skin of the satellite instead of from a retro- 
reflector assembly. However, certain important differences should be 
noted. 

First, since no retro is involved, the reference point must be 
designated by the operator during acquisition. This point should be 
selected to minimize the changes of drifting off the target, i.e. , it 
should lie within the main body of the satellite and not on one of its 
extremities. If manual tracking is employed, some experimentation may 
be necessary to determine an optimum spot size- a small spot maximizes 
the brightness and defines the best reference point, but a larger spot 
(wider beam) is less sensitive to pointing jitter. 

A second difference is that skin tracking permits a single laser 
illuminator to be used with two tracking sensors, since the spot is 
visible from different positions (this is not possible with a retro). 

On the other hand, this same property makes it necessary to use separate 
wavelengths if the two-laser rangefinder method is employed. 

Finally, reflection from the satellite skin provides a much 
weaker and more variable intensity than that possible with a retro. 
Hence, the tracker SNR will be lower and tracking or ranging may be 
degraded unless more powerful lasers are used. 

Although several implementations of skin tracking and ranging 
appear viable, a particularly attractive method is to use a single laser 
illuminator/rangefinder with a co-located rangefinder-tracker receiver 
and a second separate tracker at the opposite end of the baseline. This 
system would provide a redundant range determination by triangulation 
(two trackers) and by ranging and tracking (hybrid). 
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APPENDIX A 

MODELING A CONTINUOUS STOCHASTIC PROCESS 
BY A DISCRETE TIME STOCHASTIC PROCESS 

In this appendix, we determine precisely how a continuous random 
process is modeled by a discrete time process. The discrete time pro- 
cess can then be directly programmed on a digital computer. The primary 
interest here is in linear systems, so we restrict attention to a vector 
valued linear system. As a specific example, a scalar first-order con- 
stant coefficient linear system is modeled. 

Consider the k-dimensional continuous time vector random process 
x(t), which is modeled by [A-l]: 

x.(t) = A(t))l(t) +£(t)ujt), (A-l) 

where A(t) is a k*k matrix, B^(t) is a kxr dimensional matrix, and u(t) 
is an r-dimensional vector white noise process with zero mean and 
covariance matrix: 

E[u.{t-,) u.'(t 2 )] = Q(t-j ) $(t 2 - t-j) , (A-2) 

where 6 represents the Dirac delta function and (’) represents the 
transpose vector. 

We wish to model the continuous system by the discrete time 
linear system 

x( k+1 } = A(k)x(k) + B(k) u(k) , (A-3) 

where the dimensionality of the matrices and vectors in (A-3) is 
identical to those in (A-l). By modeling, we explicitly mean that 
the first- and second-order statistics of x(k) will be identical to 
those of xjt) at the time t=kT, where T is the sampling period of the 
discrete time process The time T can also be interpreted as the update 
time of the digital computer, whether this be a real-time or an off-line 
modeling process. 

In our model, we will make use of the state transition matrix [A-2], 
4^(t), which is the solution of the differential equation 
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l(t) = A(t)*(t) 

£(0) = £ . (A-4) 

In the special case where A(t) is a constant, i.e , A(t) = A, then 

£(t) = exp (At) , t >. 0 . 


Continuous System Model 

Equation (A-l) is a continuous linear first-order differential 
equation and has the explicit solution 


x(t) = £(t) 


)<(0) + ( £ ^ (t) £(x) U_(t) dx 
J 0 


,tiO, (A-5) 


where _x(0) is the initial value of _x{t) at t= 0. At time kT, the con- 
tinuous system in (A-5) can be written as 


x(t) = £(kT) 


x(0) + | £ ^ (t) £(x) £(x) dxj , (A-6) 


and at time (k+l)T, it becomes 


x((k+l)T) = £((k+l)T) 


x(0) + 


kT 


£**^ (t) £(x ) u_(t ) dx 


r(k+l)T , 

£ (x) £(x) £(x) dx 

kT 


. [ i r(k+l)T 

= £((k+l)T) I £ ( kT ) £( kT ) + £ (x) J3(x) u^(x) dx 


kT 


(A-7) 


Therefore, for the continuous time system, we have that 

i r(k+l)T , 

x((k+l)T) = £((k+l)T) £ _l (kT)x(kT) + £((k+l)T) £ _l (x) B.(x) u_(x) dx 


kT 


(A-8) 
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from which the expected value is given by 

E[x(k+1)T] = *((k+l)T) l _1 (kT)E [x(kT)] (A-9) 

and the covariance matrix at time (k+l)T is 

t (k+1 )T i _-] _■] 

cov[x(k+l)T] = $((k+l)T] £" (t) B_(t) (^(t) _B' (t) (-rJj'dTj*. ((k+l)T) 

J kT 

(A-10) 

Discrete System Model 

Similar expressions are obtained for the discrete time model via 
(A-3). In particular, the mean is given by 

E[x(k+1)] = A(k)E[x(k)] (A-ll) 

and the covariance matrix is 

c ov[x.(k+l)] = Bjk) E[^(k) u/ (k)] B_' (k) . (A-12) 


Equating Discrete and Continuous Time Models 

In order for the first moment of the discrete time model to be 
equal to that of the continuous model requires that E[x.(k+1)] in (A-ll) 
be equal to E[_x(k+1)T] in (A-9) for all k. This is the case when 

A(k) = £((k+l)T) ^ -1 (kT) . (A-13) 


The second-order statistics become equal when cov[x(k+l)] in (A-12) 
is equal to cov [x((k+l)T]] in (A-10) for all k. This provides the 
requirement that 


B(k)E[u(k) u'(k)]B'(k) 


= ±((k+l)T) 


f (k+l)T , 

(t) B_(t) £(t) B* (t) [ (t)]' dt ((k+l)T) 

J kT 


(A-14) 


These results are applied to specific one-dimensional cases which 
are used throughout the radar simulation. 

Examp 1 e 1 . Consider the one-dimensional version of (A-l) where 
the coefficients are constants. Therefore, 
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A(t) = A 

B(t) = B 
and 

$(t) = exp (At) (A-15) 

Also, let the correlation function of u(t) be Q(t) = (N Q /2) s(t), 
where Nq is the one-sided power spectral density of the white Gaussian 
source noise in watts/Hertz. 

Then, from (A-13), 

A(k) = exp [A(k+1)T] exp (-AkT) 


so that 


A(k) = exp (AT) , 


(A-16) 


which is a constant depending only on the parameters A in (A-14) and 
the computer update time T. 

Similarly, from (A-14), 

f (k+l )T 2 

exp (-At) B^(Nq/2) exp (-At) dx 


B 2 (k) E [u 2 (k)] = exp [2A(k+l)T] 


kT 


which reduces to 


2 2 N n B 

B^(k) E[u^(k)] = — exp (2AT) [1 - exp (-2AT)] 

4A 


(A— 1 7) 


(A-18a) 


or, equivalently. 


B 2 (k)E[u 2 (k)] 


2A 


exp (AT) sinh (AT) 


(A-18b) 


In a stable system, A in (A-14) will be less than zero. Then, 
when T + «, 

B 2 (k) E [u 2 (k)] = 0 


4 1 A| 


(A-19) 


which from (A-12) can be seen to also be the variance of x(k+l), i.e.. 
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lim E[(x(k+1)) 2 ] = N q B 2 /(4|A|) . (A-20) 

By letting T-j-® in the continuous time representation in (A-6), 
the same result is obtained. 

Example 2 . Consider the case of white noise through a first- 
order Butterworth (RC) filter, as shown in Figure A-l. 


n(t) 

w 

1 

x(t) ^ 


TS + 1 

W 


Figure A-l. White Noise Through a First-Order Filter 


The white noise has one-sided PSD of Nq watts/Hz and the filter 
time constant is x. The differential equation for the filter in 
Figure A-l is 

x ( t ) = (-1/t) x(t) + (1/t) n(t) (A-21) 

where when compared to (A-14) and (A-l), it is seen that 

A = -1/t 

B = 1/t. (A-22) 

From (A-l 6) , the constant A(k) in the discrete time model is 

A(k) = exp (-T/t) (A-23) 

and, from (A-l 8) 

9 9 M1/t) 2 

B Z ( k) E [ir(k)] = 2 (..t/- ex P (-T/t) sinh (-T/t). (A-24) 

Q 

In (A-24), either B(k) or E[u (k)J can be chosen arbitrarily. If we 
elect to set 

B(k) = 1 , (A-25) 


then 




E [u 2 (k)] 
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= ( 21 ) (t - ) exp sinh • (A-26) 

A block diagram of the model is shown in Figure A-2. The output 
of the random noise generator (RNG) is a Gaussian random variable with 
zero mean and unit variance. The sample of u(k) is obtained by scaling 
this random variable by the standard deviation of u(k). The resulting 
sequence of random variables {x(k)} will have first- and second-order 
statistics which are identical to those of x(t) in Figure A-l at times 
t = kT 

As the sampling time approaches infinity (T->-®), A{k) + 0, which 
implies that the samples are independent as is logically anticipated. 

Also, (A-l 0) shows that the variance of any sample is equal to Nq/4t. 

This is also anticipated because it represents the steady-state vari- 
ance at the output of an RC filter whose time constant is x. (The two- 
sided noise bandwidth of an RC filter B. =l/2x). 



Figure A-2. Diagram of Discrete Time Model for First-Order 

Stochastic Process 
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